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The Effects of Radiation on the near Infra-Red Absorption 
Spectrum of «-Quartz{ 
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ABSTRACT 


The variation of the absorption strength and spectrum in the region of 
3 has been determined for different specimens of crystalline quartz 
and for x-ray, electron and pile irradiated crystals. It is shown that the 
3 group of bands is associated with imperfections rather than harmonics 
or combinations of lattice vibration frequencies, and that the imperfections 
are not Frenkel or Schottky defects. Two of the bands, 2-88 and 2-95 p, 
are anisotropic, having anisotropic absorption ratios (E,,c/E,c)<0-130 and 
<0-145 respectively. A band is induced by pile irradiation at 2-79 uw but 
tends to saturate at alow dose (<10'%n°cm~?). The band is not produced 
by ionizing radiation. It is suggested that the 2-79 « band may be due to 
O-H bond formation by association of hydrogen impurity with the 
products of atomic damage. This band is also anisotropic with a ratio of 
0-71. Possible defect axis directions which could give the observed 
anisotropic ratios are discussed. 


§ 1. [INTRODUCTION 


THERE is a group of absorption bands in quartz centred about 3. The 
spectrum has been measured by, amongst others, Dreisch (1927), Plyler 
(1929), Drummond (1936) and Arnold (1954), and examination of their 
results, together with those obtained in this laboratory, indicates that 
the strength of the absorption is specimen dependent. It has to be 
assumed therefore that the bands are associated with imperfections, 
contrary to what has frequently been postulated. By introducing 
further imperfections by irradiation with high energy particles we can 
determine whether the existing bands are associated with the simple 
point defects. 
§ 2. EXPERIMENTAL 
All the measurements have been made with the cooperation of Dr. 


J. Goulden of the N.I.R.D., Shinfield, who kindly put the double-beam 
Grubb—Parson’s spectrometer at our disposal. In most cases the rocksalt 


- Communicated by Professor R. W. Ditchburn, 
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prism was used but some measurements were made using a grating 
(2400 lines per inch). A specially constructed specimen holder and cell 
enabled spectra to be obtained at temperatures in the range —190°c to 
100°c. 

Most of the measurements have been made on specimens of the best 
optical quality Brazilian quartz obtained from Hilger and Watts, Ltd. 
The specimens were 1:5X1:01:0 cm and had two pairs of polished 
faces so that anisotropic absorption measurements could be made on the 
same crystal. The c-axis was along one of the 1cm edges. Some 
measurements were made on a synthetic quartz crystal obtained from 
the G.E.C. Research Laboratories, Wembley. 

The absorption of | cm of quartz at 2 p is negligible (less than 0-01 cm, 
Drummond 1936) and from the known refractive index the transmission 
of a non-absorbing specimen for light at normal incidence is calculated 
to be 91:5°%. Between 2 and 4y the refractive index can be taken as 
constant without introducing any significant error. The double beam 
compensation mechanism was correct to -+-0-5°% for transmissions 
between 10 and 100% so that the value of transmission (7’) at a given 
wavelength was found by linear interpolation between the zero and 
91-59% levels. If «cm is the specimen thickness the absorption co- 
efficient is given by u(cm-!)=1/z In 91-5/T. The errors in » resulting 
from the measurement of the traces are --0-1 cm~! for a 1 em specimen 
when »=2-5 cm™!, falling rapidly to +0-01 em—! for p.20-1 em-!. The 
wavelength scale was calibrated using the lines of water vapour and 
from a polystyrene sample using the wavelengths given by Plyler and 
Peters (1951). 

Details of irradiations are given in the appropriate sections below. 


§ 3. RESULTS 


3.1. The Absorption in Unirradiated Crystals 


The absorption spectra of various unirradiated crystals are shown in 
figs. | and 2. All the curves refer to measurements in which the light 
travels parallel to the c-axis. It will be seen that these spectra are com- 
posed of several overlapping bands which we describe as «, B, y ete. 
Labelling the bands in this way is convenient because the value of the 
photon energy at the peak depends on the amount of absorption present 
in the adjacent bands. 

In some specimens (fig. 1) the «-band is prominent, whilst in others 
(fig. 2 (a)) it is barely resolved. The band £ is just a point of inflection 
on all the curves but y, 6 and ¢ are well defined peaks. The y-band is 
the principal peak in all but one of the specimens and occurs at a photon 
energy of 0-420ev (2-95). The exception is the synthetic crystal 
(fig. 2 (b)) in which both « and y are weaker, possibly absent, and a peak 
which may correspond to f, is resolved between the normal «, y positions, 
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Absorption curves of unirradiated natural crystals («-band strong). 
A, Mitchell and Rigden, B, Drummond, C, Dreisch. 
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Absorption curves of unirradiated natural and synthetic crystals : 
(a) Natural crystals, «-band weak ; ——-—————Mitchell and Rigden, 
—, — — —Arnold. 
(b) Synthetic crystals: curves A and B, Mitchell and Rigden A-E, 
c-axis, B-E, c-axiS; curves C and D, Arnold. (Assuming Arnold’s 
ordinate to be in cm-!.) 


944 E. W. J. Mitchell and J. D. Rigden on the Effects of 


Synthetic quartz sometimes has a broad absorption band at 3-ly 
which has been attributed by Brown et al. (1952) to the presence of water 
inclusions. They found that crystals having this band also showed 
considerable Tyndall scattering of visible light, the scattering presumed 
to come from the inclusions. The synthetic quartz measured in the 
present work was cut from one of the high quality, low aluminium con- 
centration (ca. 0-01 atom%,) crystals now grown by Brown and Thomas 
(1956) and does not show the water inclusion band. Neither has this 
band been found in any of the natural crystals which we have examined. 


Fig. 3 
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Although free from water inclusions the synthetic specimen (fig. 2 (b)) 
did have a peak at 2-81 4 which was not present in all the natural 
specimens. It seems likely that this corresponds to vibrational absorp- 
tion of isolated O—H bonds as will be discussed below. 

Increasing the resolving power of the instrument by using the gratin 
improved the resolution of the peaks slightly and revealed some fine 
structure. Even using the rocksalt prism, however, the bands sharpened 
considerably when the specimen was cooled to —193°c as shown in fig. 3 
At the low temperature f is clearly resolved and « appears as a ae 
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3.2. Anisotropy of Absorption 

Drummond (1936) reported that the absorption around 3 ~ was aniso- 
tropic but his measurements were not made on the same specimen. The 
anisotropy is clearly shown in fig. 4, in which curve (A) refers to a measure- 
ment with light travelling parallel to the c-axis, that is with E always 
perpendicular to c, and curve (B) to light travelling perpendicular to the 
axis so that there are components of E parallel to c. 

It will be seen from fig. 4 that Me > He, where « is the absorption 
coefficient. The ratio of these quantities can be determined without 
using polarizing devices if it can be once established that the absorption 
is isotropic about the c-axis and that the wavelength at the peak does 
not depend on the direction of E. 
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The anisotropy of the absorption bands in unirradiated natural quartz. 


The absorption for all directions of E perpendicular to ¢ was determined 
as follows. Light falling on the specimen was plane polarized by a 
5 plate selenium polarizer. The crystal was mounted in a cell which 
could be rotated by a small electric motor, care being taken to ensure 
that the aperture was constant throughout the rotation. For wave- 
lengths of 2-0 and 2-95 the transmission was recorded as the specimen 
rotated. By rotating the crystal rather than the polarizer it was not 
necessary to make any correction for-the variation of efficiency of 
the monochromator and detector for different directions of E. This 
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experiment showed that the transmission was constant to better than 
4° for all E perpendicular to c. It can be assumed therefore that the 
transitions are dipolar. That the position of the peak is independent 
of (E) can be shown to sufficient accuracy by measuring with unpolarized 
light passing along and perpendicular to c. 

We can now determine the transmission of the crystal for E,c (7',) 
and for E,c (7) without using a polarizer, and independently of any 
selective polarization effects in the monochromator or the detector. 
Let [,(O) be the intensity of unpolarized light incident on the crystal 
in the direction E,c, J, the intensity transmitted by the crystal assuming 
that the reflection correction has been made, J ,(O) and J, the correspond- 
ing intensities when the light travels in a direction (1) perpendicular to c, 
I,(O) and J,, for a direction (2) perpendicular to ¢ and (1), and k, and ky 
be the sensitivities, relative to some arbitrary direction, of the mono- 
chromator and detector for E parallel to (1) and (2) respectively. Then 


a key + al kes f, 


tee 47,(O)k,+41,(0)k, I,(0) i 
U eee 41 (O/T kh, +41,(O)T ke = kT +kT\ 
2 41,(O)(ky +k) ki +k, 


koT'+k,Ty 


M— gM eee ot Ae) = (37,+T\) 
and T \=22 4-7 
For the most accurate measurements it is essential to use a beam of 
parallel light whereas in our experiments the beam was diverging. The 
main uncertainty in the measurements to be described, however, does 
not arise from the divergence but from the presence of overlapping bands. 

Using this technique curve (C) in fig. 4 was obtained for the absorption 
spectrum for E,c, the perpendicular case still being given by curve (A). 
The difference curve (D) shows that both «- and y-bands are strongly 
anisotropic and that 6 is probably anisotropic. The « and y peaks are 
still present in the E,c curves but owing to the overlap of the bands 
and the uncertainty about the 8-band we can only put a lower limit to 
the ratio (.) of the absorption when the E is parallel to c, to that when 
E is perpendicular. Taking the peak values we find .,—0-130 and 
@,=0-145. If it is assumed that the absorption for E,c in the region 
of the « and y peaks is composed only of « absorption and y absorption, 
then resolving geometrically we find .7,—0-110 and &~,=0-123. 

The final curve (EK) in fig. 4 was obtained using the selenium polarizer. 
Using the polarizer alone the deflection from the detector was greater 
when E was horizontal than when vertical. Accordingly the polarizer 
was mounted to give E horizontal and the specimen mounted so that the 
light passed perpendicular to the axis with E,c. The difference between 
curves (EK) and (C) is due to the polarizer not being 100°% efficient, the 
effect of the unpolarized component being enhanced by the high aniso- 
tropic ratio. 
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3.3. The Effects of Irradiation 

Although the study of the effects of irradiation is made complicated 
by the overlapping band structure there are two unequivocal results 

(a) None of the bands «~e is continuously enhanced by pile irradiation. 

(b) A new band (x), which was briefly reported by Paige (1955), is 
induced at 0-444 ev (2-79,) by pile irradiation, the value of the peak 
energy at any time depending on the amount of overlapping from the 
lower energy bands. 
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The effect of pile irradiation on the near infra-red absorption of natural quartz : 
curve A, unirradiated; B, neutron dose 6-510" 7° cm-*?; and C, 


neutron dose 44 « 1017 n° em-?. 


These two results are illustrated in fig. 5 in which absorption curves 
are plotted for two neutron irradiations. The visible colouring of this 
specimen due to aluminium impurity, Mitchell and Paige (19955), was 
at all times uniform. 

It is also clear from fig. 5 that both « and y decrease rapidly with 
pile irradiation. The extent to which these changes are due to ionization 
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effects has been investigated using ultra-violet illumination, x-ray and 
2 mev electron irradiation. Illuminating an untreated crystal with ultra- 
violet light from a high pressure mercury lamp had no effect on the 
absorption spectrum. Two opposite polished faces of another crystal 
were irradiated with x-rays from a copper target, the dose being ca. 1071 
per face. This treatment had the effect of decreasing the «-band such 
that it was no longer resolved. The electron irradiation was carried 
out at the Tube Investments Research Laboratory, nr. Cambridge, using 
the 2mMev Van de Graaff accelerator. The macroscopic temperature 
of the specimen during the irradiation was —80°c, the total electron 
dose being 950 (uA mincm-2) or 3:6 1017e cm-?. This produced 
a very intense visible absorption band (A,, A,—associated with Al, 
Mitchell and Paige 1955) in the crystal and led to the decrease of absorp- 
tion in the we bands. The effects of ionizing radiation on the 3 « bands 
is exceedingly complex and further work is being carried out. 

Approximate calculations, and the evidence of the visible colouring, 
show that as far as ionizing effects are concerned the above electron 
irradiation would give more ionization in the crystal than would have 
occurred during the pile irradiation of these specimens. It is concluded, 
therefore, that the y-band, which is induced only by pile irradiation, is 
associated with atomic displacements. 

The neutron dose corresponding to curve (B) in fig. 5 would be expected 
to induce about 10!%cm-* displacements, Mitchell and Paige (1956), 
whereas for a displacement energy of 40 ev the 2 Mev electron irradiation 
would be expected to give of the order of 101*6cm-%. Thus, although 
the y-band is associated with displaced atoms, because its strength is 
small after pile irradiations we should not expect to see it after the electron 
irradiation. 

Since the y-band is associated with an atomic displacement the rate 
curve for the band is of considerable interest. For the neutron doses 
used during the present experiments (up to 50x 1017 n° em-2) the rate 
of replacement collisions can be neglected. For effects associated with 
simple point defects we should therefore expect an approximately linear 
increase with neutron dose. The measured rate curve is shown in fig. 6 
and is clearly not linear, the slope continuously decreasing. In addition 
measurements have been made using a stack of five 0-2 mm. thick speci- 
mens which had received an estimated neutron dose of 80 1017 n° em-?. 
The y-band could just be detected. These specimens were irradiated 
for another purpose and were too small for accurate infra-red measure- 
ments. However, they show that the y-band absorption is less than 
1 cm~! for this dose. 

Because of the decrease of the «- and y-bands the real rate of increase 
of the y-band near the origin cannot be determined exactly. The lowest 
rate (A) is given by joining up the experimental points for the peak 
absorption and an estimated real rate, allowing for the decrease of the 
a- and y-bands, is also shown (B). 
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In order to try to decide the axis of the y centre we have determined 
the anisotropy of the y-band. This was not possible with any accuracy 
until after the dose of 44 10!7 7° cm-2. Then the «- and y-bands were 
sufficiently small and the radiation induced band (x) sufficiently large 
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Rate curve for the y-band induced by neutron irradiation. (Curve B is 
corrected for the decrease of «, assumed to be of the form C). 


for a specific value of 7, to be measured. The curves are shown in 
fig. 7, from which a value of ,=0-71 is obtained. The same reasons 
which made these curves suitable for determining /, make them un- 
suitable for determining further values of 7, and .,,. 


§ 4. Discussion OF RESULTS 


4.1. The Origin of the 3. Group 

Several authors have suggested that these bands are harmonics and 
combinations of the fundamental lattice vibration frequencies. By 
choosing suitable combinations Plyler (1929) predicted wavelengths at 
which such bands would be expected to occur and some of the predictions 
corresponded to the observed peaks. Saskena (1940) in his detailed 
analysis of the vibrational spectrum of quartz comments that the 3 yu 
bands probably result from harmonics and combinations of the lattice 
vibrations. Arnold (1954) made the specific suggestion that the y-band, 
which he found at 2-92 , is the third harmonic of the 8-7 . Si-O vibration. 

This type of explanation cannot be maintained, however, because 
of the wide variations in absorption strength of the bands in different 
crystals. In order to account for these variations it would be necessary 
to assume that about one half of the bonds were broken in some specimens 
during growth of the crystal, and in others, as the result of relatively 
weak irradiations during which only 10'%-101!9 cm~? defects would have 
been produced. Neither assumption is valid. In any case Arnold’s 
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suggestion is ruled out by the large anisotropy of the y-band 
(He, Cm™*>/4g,) compared with yg, =0 for the 8-7 » band. 

It is necessary to assume, therefore, that the bands are associated 
with imperfections. None of the bands is progressively increased by 
pile irradiation—cf. the C-band, Mitchell and Paige (1956)—indicating 
that the original bands are not associated with displaced silicon or oxygen 
atoms, or the vacancies. Consequently it is likely that the defects 
responsible for the bands are impurity atoms, or are in some way associated 
with the presence of dislocations. 
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Anisotropy of the y-band. 


The measurements at low temperature showed that both « and y 
decreased in total strength as the temperature was lowered. For electro- 
nic transitions involving at least one localized level in a solid the total 
absorption strength is generally independent of temperature. The 
observed decrease can be qualitatively accounted for in several ways: 

(a) A normally forbidden electronic transition of an impurity, the 
atomic selection rule being relaxed in the crystal as a result of lattice 
vibrations—e.g. Cu in NaCl, Seitz (1951). 

(6) A centre having three electronic levels, L, M and N such that the 
L-M transition is allowed but L—N forbidden. If there are in the ground 
state two electrons in each of L and M no absorption will be possible. 
When electrons from M have been thermally excited to N the L-M 
transition will be observed—the higher the ‘population of N (higher 
temperature) the higher the strength of L-M absorption. ‘ 

(c) A vibrational absorption associated with an impurity. 
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Regarding (c), the energies involved (ca. 0-4 ev) are high for impurity 
induced vibrations, and the absorption strength is too high for the 
bands to correspond to harmonics of a lower energy impurity induced 
vibration. 

The normally forbidden transition (@) becomes allowed because in the 
real crystal the atomic wave functions (e.g. S,, S,) are modified, having 
to a limited extent the character which would give an allowed transition 
(e.g. S,, (S+dp),). It is this which gives the transition probability of 
the observed band and we should expect the anisotropy to be character- 
istic of dipole rather than quadrupole transitions. An explanation of 
type (a) cannot be ruled out, therefore, on the grounds that « and y 
absorption is isotropic for all directions of E perpendicular to the c-axis. 

The centres giving rise to the 3. bands must have the property that, 
on capturing or losing an electron, the absorption—vibrational or electro- 
nic—is drastically reduced. For a centre of type (b), for example, a 
trapped electron would either go into N thereby inhibiting absorption 
or the levels would be modified giving absorption at some longer wave- 
length (cf. F-F’ in NaCl), with a consequent decrease at 3p. 

To account for the differential changes in absorption which have 
been produced by cooling, by using polarized light and by irradiating with 
ionizing radiation it is necessary to assume that more than one type of 
centre must be contributing to the 3 absorption bands. 


4.2. The y-Band 


It has been shown that the radiation induced band is associated with 
atomic displacements. The departure of the rate curve from linearity 
indicates that the rate of displacement of atoms cannot be the complete 
rate determining process. We conclude from this and from the tendency 
to a saturation value of absorption (fig. 6) that a product of damage is 
becoming associated with an impurity (or dislocation). 

The peak absorption occurs at 2-79 (--0-01) which is close to the 
wavelength at which the O-H bond vibrational absorption occurs in 
vapours and weak solutions—e.g. 2-696 to 2-770 u, Fox and Martin (1937), 
(1940). A strong band has been found in some specimens of fused 
quartz, and in glasses, at 2-75 and this band has been attributed to 
O-H vibrations—e.g. Ellis and Lyon (1937), Florence e¢ al. (1953). 

Quartz probably crystallized in nature by a hydrothermal process 
similar to that now used for producing synthetic quartz. Consequently 
hydrogen is a likely impurity and O-H bonds could be formed during 
irradiation by the combination of interstitial hydrogen and interstitial 
oxygen, or by hydrogen occupying silicon vacancies. The concentration 
of O-H bonds which would be required to account for the observed 
absorption strength can be calculated using the probability of absorption 
determined from the measurements on weak solutions. Kaye (1954) has 
given a mean value of the peak absorption per O-H bond in various solu- 
tions, and assuming an average width at half-absorption of 0-04 ev—based 
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on the measurements by Fox and Martin (1937), (1940)—we find a 
mean value of (2-+-1) x 1021 O-H bonds (cm) per unit integrated absorp- 
tion expressed as (f,,,< 4E) in cm~1 ev. The initial rate of increase 
of the y-band, allowing for the effect of the decrease in the «- and y-bands 
on both the rate and the bandwidth, is 2x 10-3 cm~! ev per 1017 n° cm-?. 
This would be equivalent to the formation of 4 1018 (em~%) O-H bonds 
per 1017 incident neutrons per cm?. Mitchell and Paige (1956) estimated 
that ca. 3x 1018 displaced oxygen atoms and 1x 1018 displaced silicon 
atoms are produced per 1017° cm~*. If the rate of formation of O-H 
bonds was limited further by diffusion of the hydrogen or the displaced 
atoms then the observed rate should correspond to a lower limit for the 
rate of production of point defects. To an order of magnitude, therefore, 
the above interpretation is consistent with the observations. Further 
experiments are being carried out to determine whether the hydrogen 
concentration in quartz can be varied in a controlled manner by heat 
treatment, thus enabling more quantitative tests of the mechanism to be 
made. 

Whatever the origin of the y-band, the band which was resolved at 
2-81+0-01 in unirradiated synthetic quartz (curve 1, fig. 2(c)) is 
likely to be associated with O-H. This band is also resolved on one of 
Arnold’s curves for synthetic quartz, while in unirradiated natural 
crystals we have sometimes detected changes of slope in absorption 
curves which could be associated with a band at about 2:80. However, 
the anisotropy differs from that of the y-band. For the y-band 
pcm *y,,em~1 is 1-4 while for the 2-81 1 band in unirradiated synthetic 
quartz j,/4, is at least greater than 2-2 (see fig. 1 (c)). It is possible, 
owing to the way in which the crystal was grown artificially with the 
major growth perpendicular to the c-axis, that O-H groups tend to be 
incorporated preferentially in the growth plane. 


4.3. The Anisotropy of the «-, y- and x-Bands 


For a set of 3n linear oscillators regularly arranged about the c-axis 
with which they make an angle 6, the anisotropic ratio ./ is 2 cot? @. 


If we have more than one set 371, 3n....... making angles 0,,0,..... 
then 
2 
Per >,3n, cos? Ox 


fa”, 8/2 Sino 03 Sewn) 


Thus for defect oscillators having axes along Si-O directions for which 
G,=44° 18’, 0,=66° 15’, ny=n,=2 the ratio ~W is 1:06. For the other 
simple directions, parallel to the following, the values are: c, =o; 
a, @—0; nearest neighbour Si-Si, .7=1-06; perpendicular to Si-Si 
through the common O, ..=0-74 ; and perpendicular to the plane of the 
Si-O-Si triangles, .7=1-3. 

The anisotropic ratios which have been measured in the infra-red 
are: %, 0-130; y, 0-145 and y, 0-71. If each of these is interpreted as 
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being due to sets of oscillators at one angle 0, we find 075° 42’ for the 
a-band, 74° 57’ for y and 59° 12’ for the y-band., 

As far as « and y are concerned there is the possibility that the oscillators 
lie in planes parallel to (001)—e.g. linear oscillators along the a-axis 
tunnels—and that the measured anisotropy differs from zero because of a 
combination of beam divergence, misorientation of the crystal and 
thermal vibration of the axis out of the plane. 

The beam of light used in the experiments was diverging with a maxi- 
mum semi-angle of 6=2°. Thus for light nominally travelling perpendi- 
cular to ¢ and having E nominally parallel to ¢ we have contributions 
Esin¢ capable of exciting the centres. In the case of E nominally 
perpendicular to c the effective E becomes E cos ¢, so that for a prism 
shaped light path in the crystal 

2 

t= | ik tan? ¢ dd, 

=o 
which is negligible for a divergence of +2°. Possible errors due to 
misorientation, including cutting and mounting, were less than 2°. 
The effect of temperature on the anisotropy cannot be easily determined. 
The difficulty is that the uncertainty in the anisotropic ratios, leading 
to the values of 0-123 and 0-145 being considered as upper limits, arises 
_ because of the uncertainty in determining the amount of « and y present 
in the E,c case when both « and y are small and overlapping. At low 
temperatures the « and y bands are reduced and this makes the determina- 
tion of the anisotropic ratio more difficult. 

The ‘ thermal spread ’ suggestion cannot be ruled out experimentally 
at present, although the angular variations required to account for the 
results (a mean of +15° to give observed upper limit of .—0-13) would 
have to correspond to approximately 0-5 A displacement out of the plane. 

There are two directions between well defined points in the structure 
which make about the required angle to the c-axis. Oxygen atoms (9), 
(5), (4w) and (7b)—using the numbering of atoms shown on fig. 8—form a 
tetrahedron, and a similar tetrahedron is formed by the oxygen atoms 
(7b), (Sb), (6b) and (4). If P and P’ are respectively the points which 
are equidistant from the atoms in each of these tetrahedra, then PP’ 
makes an angle of 67° with the c-axis. For an oscillator in such a direction 
f =0-36 compared with the observed upper limits of 0-130 and 0-145, 

The second direction is one of the O—O directions. The O—O directions 
are: 6 similar to (5) to (8), 78° 50’ with c, oxygens separated by 2-60 A; 
6 similar to (5) to (6), 47° 3’ with c, 2-64 A ; 3 similar to (6) to (8), 60° 44’, 
2-65 A; and 3 similar to (5) to (7), 27° 49’, 2:60 A. Taken together and 
using eqn. (1) the anisotropic ratio is 1:02. Taking the directions between 
the closest oxygens—viz. directions making 78° 50’ and 27° 49’ with c¢, 
the ratio is 0-80. By itself the 78° 50’ O-O direction gives a value for 
the anisotropy of 0-078, which is consistent with the observed limits. 
It is not clear, however, if the centre is diatomic, how the axis is parallel to 
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this particular O-O direction, or if the centre is a simple impurity how 
the directed wave function has its axis in this direction, rather than being 
distributed about the four O-O directions. The same difficulty was 
noted by O’Brien (1955) when explaining the anisotropy of some of the 
paramagnetic resonance effects associated with aluminium impurity in 
quartz. To account for the quadrupole effects Miss O’Brien had to 
assume that only two of the possible four oxygen atoms surrounding a 
given aluminium were suitable for the formation of the colour centre. 


Projection of atom positions in « quartz on to a (001) plane. Atoms in the 
unit cells below those shown are given a suffix b. 


The y-band appears to be related to the association of a product of 
damage with an impurity. The simplest possibilities are shown in the 
table. We have omitted from this table the combinations of an impurity 
substituted for silicon with interstitial silicon or oxygen atoms, and of an 
impurity substituted for oxygen with interstitial silicon or oxygen atoms. 
There are large numbers of such directions and a random distribution of 
defect-pairs will tend to be isotropic. 

The values of anisotropic ratio given in the table were calculated 
from eqn. (1). In calculating the ratio for interstitial pairs we have 
assumed that these are distributed between the major tunnels in the 
structure—those parallel to the c-axis; and those parallel to a-axes. 
For the rhombohedral unit cell there are 3 a-axis tunnels per cell, and 
1 c-axis tunnel. The c-axis tunnels are 1-1 times longer than those 
along the a-axes so that we take .7=1-1/(3/2)=0-73, in good agreement 
with the measured value of 0-71--0-04. All the other values given in 
the table differ significantly from the observed ratio. 

There is, however, a further possibility to be considered that the 
observed transition is associated with a wave function directed at right 
angles to the principal axis of the centre. For example, the tightly 
bound electrons of the centre may have p-function components parallel 
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to the principal axis, while the observed transition may be associated 
with a less tightly bound electron having a p-function perpendicular 
to the main axis. In particular, a centre comprising a substitutional 
impurity in a silicon position adjacent to a silicon vacancy, would give an 
anisotropic ratio of 0-74 for the perpendicular direction passing through 
the nearest oxygen. This direction can also account for the observed 
ratio. 


Anisotropic ratios for some impurity/damage combination centres. 
Where the centre can occur equally in several directions values 
are given for the mean and also for each of the directions 
assuming that all the centres are in one direction. 


Position of Damage Principal Anisotropic 
impurity product axis, 0 ratio 
Substituted for Si Si vacancy 54° 1-06 
O vacancy 44° 18’ (2-10) 
1-06 
66° 15’ (0-39) 
Substituted for O Si vacancy 44° 18’ (2-10) 
1-06 
66° 15’ (0-39) 
O vacancy 27° 49’ (7-18) 
ime (1:75) 
1-02 
60° 44’ (0-63) 
78° 50’ (0-08) 
Interstitial Interstitial 0° 0:73 (00) 
Si or O 90° (0) 


§ 5. CONCLUSIONS 


The experiments have shown that the 3 1 group of bands is associated 
with imperfections in the crystal and that the imperfections concerned 
are not interstitial silicon or oxygen atoms, nor silicon or oxygen vacancies. 
At least two centres are required to account for the various differential 
changes in absorption strength. The absorption in the 3 u bands is 
modified by the electron and hole trapping which accompanies the 
passage of ionizing radiation through the crystals. 

Two of the bands are very anisotropic and the important axes associated 
with the observed transitions are 75° 42’ («) and 74° 57’ (y) from the 
c-axis, assuming that all the oscillators of each kind make the same angle 
with the c-axis. This could correspond to a direction parallel to one of 
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the nearest neighbour O-O directions but there is no other evidence 
to show how one such direction should be preferred to the other three 
O-O directions. 

The radiation induced band is associated with atomic displacements 
and their combination with some imperfection—possibly hydrogen 
impurity. The anisotropy of this band can be accounted for if the 
centre is a combination of an interstitial oxygen, or silicon, atom with 
an interstitial impurity. Thus the anisotropy is consistent with the 
O-H model. Alternatively, agreement with the observed anisotropic 
ratio is found for an oscillator having an axis in a direction perpendicular 
to a line joining two neighbouring silicon positions and passing through 
the nearest oxygen atom. 
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ABSTRACT 


A xenon filled proportional counter has proved to be a very useful 
instrument for the study of z-mesonic x-rays from light elements. 
Measurements of the K, energies from *Li, 7Li, Be, B and “B are 
described. Energy shifts due to the short range pion-nucleon forces 
are observed, confirming the scintillation counter measurements. 

The experimental values of the energy shift are fitted equally well by 
the theories of Deser et al. (1954) or Karplus (1956). 

Measurements of L, energies are used to deduce a value for the 7-meson 
mass equal to 272-5-+-0-8 m,. 

The natural broadening of the beryllium K, line due to the short lifetime 
for capture of the 7-meson from the Is state has been measured. A width 
['=1-2+0-2 kev was observed. 

Intensity ratios of K, and Kg and L,, Lg and L, lines are also reported. 


§ 1. INTRODUCTION 


7™-MESONIC X-rays were first detected by Camac, McGuire, Platt and 
Schulte (1952). Later work by McGuire et al. (1954), and by Stearns, 
DeBenedetti, Stearns and Leipuner (1954), established that the yield 
of the K series x-rays decreased rapidly with increasing Z due to capture 
from the 2p level, and was barely detectable above Z=9. 

The first accurate measurements of the energies of the K x-rays were 
made by means of critical absorbers in conjunction with a measurement of 
the pulse height distribution from the x-rays with a sodium iodide 
scintillator. Stearns, Stearns, DeBenedetti and Leipuner (1954), showed 
that the energies of the K, X-rays from lithium and_ beryllium 
were definitely lower than the values calculated from purely electro- 
magnetic interactions. Later, improvements in the technique enabled 
direct measurements of energy to be made from pulse height distributions 
obtained with a Nal scintillator. Stearns, Stearns, DeBenedetti and 
Leipuner (1955), Camac, McGuire, Platt and Schulte (1955) and Stearns 
and Stearns (1956 a) measured the energies of 7-mesonic K x-rays for 
elements up to fluorine and obtained energy values which were in all 
cases lower than those calculated from purely electromagnetic interactions. 
The energy shifts are due to the short range forces between pions and 
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nucleons and the observed shifts were in qualitative agreement with a 
theory of Deser et al. (1954) discussed below. 

In the case of the lighter elements the energy resolution of the scintilla- 
tion counter is not small enough to resolve the K , (2p— 1s) and Kz (3p 1s) 
lines. This is well within the capabilities of the proportional counter 
whose resolution is at least four times narrower than that of the scintilla- 
tion counter. It was therefore desirable to repeat these measurements 
with a proportional counter. Moreover the improved resolution enabled 
the natural broadening of the 1s level in beryllium to be measured for the 
first time. The level broadening arises from the short lifetime for 
capture of the 7-meson in the 1s state. Measurements of the energies 
of L, x-rays were also carried out and the relative intensities of K, and 
Kg and also L,, Lg and L, radiations can be estimated. 

The first results obtained with proportional counters were described 
in a letter, West and Bradley (1956 a) and interim accounts of the present 
work were presented at the 6th Rochester Conference 1956 and the 
CERN Symposium 1956 . 


§ 2. THE PROPORTIONAL COUNTER 


The two main requirements of the proportional counter were a high 
efficiency for detecting x-rays combined with a low background counting 
rate. It was therefore made capable of withstanding pressure and was 
filled to 3 or 5 atmospheres pressure with xenon. To reduce the back- 
ground counting rate due to general radiation from the cyclotron, the 
active length of the counter was kept to a minimum, which required 
the use of field tubes of the type described by Cockroft and Currant (1951). 
A diagram of the counter is shown in fig. 1. It was made of aluminium 


Fig. 1 
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Diagram of the proportional counter. 


and opposite the active region of the counter the wall was turned down 
to 0-030 in, to minimize absorption of x-rays. The field and guard 
tube assemblies were of the type described by Batchelor et al. (1955). 
They were bonded with ‘ araldite’ to invar plugs which were screwed 
and soft soldered into the steel core of a section of ‘ Al-fin ’ the aluminium 
outer portion of which was argon arc welded to the counter body. The 
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inside of the counter wall was finely polished and then oxidized in the 
manner described by Wilkinson (1950) for Geiger counters. 

It had previously been noticed that the uniformity of the counter wire 
was very important in pressurized proportional counters. The wire 
used here (0-005 in. Cu-Ni alloy) was specially prepared by Johnson 
and Matthey Ltd. using a freshly prepared diamond die. The uniformity 
of multiplication factor was tested by measuring the pulse height distri- 
bution due to the 59-6 kev y-rays from *41Am with the source placed at 
different positions along the counter and round its axis. With a filling 
consisting of 3 atmospheres pressure of xenon +3 cm of Hg pressure 
of methane, the counter gave a multiplication factor which was uniform 
over the active volume to better than half of one per cent when the field 
tube potential was correctly adjusted. 

The counter operated with an H.T. voltage of ~ 5000 volts on the 
wire and amplifier integration and differentiation time constants of 
34.sec were used. It was necessary to embed the external portion of 
the field tube assembly in ceresin wax to eliminate electrical breakdown. 
The energy resolution, measured with 59:6 kev y-rays was 2:9% (half 
width at half maximum) and is only slightly worse than the best resolu- 
tion obtainable (2-5°%,) with low pressure fillings. The filling gas was 
occasionally repurified over heated calcium, or over a sodium mirror, 
during the series of measurements when the energy resolution had 
deteriorated slightly. 

Fig. 2 
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Diagram of the experimental arrangement. 


§ 3. EXPERIMENTAL METHOD 

Negative mesons from the Liverpool Synchrocyclotron entered the 
apparatus shown in fig. 2. Counters 1, 2 and 3 are plastic scintillators 
2-5 in. x 2-5 in. X 0-25 in. viewed by 2in. photomultipliers in contact 
with a small face of the scintillator. Counter 4 is a plastic scintillator 
4in.X3-5in. 0-15 in. It was used as an anticoincidence counter and 
to obtain a maximum efficiency it was viewed by two photomultipliers 
in contact with the small faces at opposite sides of the crystal. The 
two multipliers had independent H.T. supplies and their outputs were 
mixed via cathode followers. The anticoincidence scintillation counter 
completely covered the active region of the proportional counter. 
7-mesons which stopped in the target material were selected by means of a 
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coincidence between the output of counters 1, 2 and 3 (resolving time 
0-3 psec) in anticoincidence with the output of counter 4. (123-4 event.) 

The fraction of these events which gave a pulse in the proportional 
counter was further selected by means of a triple coincidence between 
counters 1, 2 and 3 in coincidence with a proportional counter pulse and 
in anticoincidence with the output of counter 4. (123 Prop.—4, event.) 
An event of this type opened a gate 11 psec in duration which allowed 
the proportional counter pulse to pass to a 70 channel pulse analyser 
of the Hutchinson Scarrott type. 

Tests with «, y coincidences from *41Am, detecting the «-particles with 
a scintillation counter and the y-rays with a proportional counter, showed 
that a maximum coincidence rate was obtained with 2 psec resolving time 
when the «-particle pulses were delayed by 2psec. Accordingly the 
triple coincidence pulse from the scintillation counter telescope was 
delayed by this amount before mixing with the proportional counter 
pulse. (With a 5 atmosphere filling 3 usec delay and 4 psec resolving 
time were used.) 

The negative meson beam consisted of 7-mesons of energy 96 Mev 
together with 6% of y’s and 4°% electrons having the same momentum 
and therefore a longer range than the 7’s. An absorber was placed in 
the beam and its thickness adjusted to give a maximum number of 
m-mesons stopped in the target. Polythene was used to minimize 
scattering and the last few inches were inserted between counters 2 and 
3 enabling the 12 coincidence rate to be used as a monitor when deter- 
mining the stopped meson rates. The mesonic X-ray sources were 
23 in. x 22 in. with a thickness of ~2 g/cm? which was known to give 
only a small amount of backscattered Compton radiation. When 
powders were used they were compressed into thin aluminium boxes 
weighing a few grammes by means of a hydraulic press. With the 
optimum absorber thickness, approximately 0-6 of the 7-mesons were 
stopped by the target material. 

In addition to 12 feet of concrete and steel shielding between the 
cyclotron and experimental rooms, a cave of concrete and iron 2 ft 
thick and open at the rear was constructed round the apparatus. The 
meson beam was magnetically deflected so that direct radiation from the 
cyclotron could not pass through the apparatus. 

The counting rate which could be accepted was limited by the single 
rate in the proportional counter. Above about 6000 c¢.p.m., which 
corresponds to a mean rate of one proportional counter pulse per cyclotron 
pulse, some broadening and additional background was observed in the 
spectra. This was presumably caused by pile up of consecutive pulses 
and throughout the measurements the proportional counter single rate 
was never allowed to exceed 6000¢.p.m. The duration of the cyclotron 
pulse was ~ 100 sec and as an additional precaution a paralysis time 
of 100 «sec was put on the gating pulse and in the later measurements a 


paralysis of 200sec was put on the proportional counter channel of 
the coincidence unit. 
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Table 1 gives a typical set of counting rates used in a measurement of 
the beryllium K, line together with the corresponding values in a 
scintillation counter experiment (Stearns and Stearns, private com- 
munication). 

The overall efficiencies of the proportional counter and scintillation 
counter are comparable, but the counting rates which can be accepted 
are much smaller in the case of the proportional counter. The back- 
ground of very large pulses in the detectors is much smaller with the 
proportional counter presumably due to its smaller stopping power. 
A major part of the background was associated with the (123-4) rate 
obtained without the target. This rate was usually about 0-5 to 0-75 
times that due to the target. It arose from mesons which stopped in the 
third scintillator and in the part of its photomultiplier and case which 
was exposed to the beam. 


Table 1. Comparison of Proportional Counter and Scintillation Counter 
Experiments 


Proportional counter | Scintillation counter 


Coincidence resolving time 


in the detector channel 2 usec 0-05 usec 
Gate length 11 psec 2:0 psec 
Detector single rate 4100 ¢.p.m. ~ 12000 ¢.p.m. 
Stopped meson rate 1070 ¢.p.m. 16000 ¢.p.m. 
Gating rate 41 ¢.p.m. 900 ¢.p.m. 
Estimated casual rate 2-6 ¢.p.m. <20 ¢.p.m. 
Gated pulses comprise : 

Mesonic X-rays 16 c.p.m. ~ 160 c.p.m. 

Background on scale 12 ¢.p.m. ~ 140 c.p.m. 

Background off scale 13 c.p.m. ~ 600 ¢.p.m. 


Energy calibrations were carried out with radioactive K x-ray sources. 
The proportional counter pulse was made to produce a gating pulse by 
feeding it into the channel usually fed by the triple pulse from the counter 
telescope as well as into its own channel, the anticoincidence counter 
being disconnected. The source was placed midway along the counter, 
the uniformity measurements already referred to having indicated that 
the multiplication factor for a source at this point differed by only a 
very small amount, <0-1°%, from the mean value over the whole effective 
region of the counter. 


§ 4. ENERGIES OF 7-MEsonic K, X-Rays 
Runs of up to 10 hours duration were needed to obtain an adequate 
number of counts in the spectra of K, x-rays. As a check on drifts 
in gain during these long periods, energy calibrations were carried out 
every two hours. The standard deviation of the mean value of the 
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position of the calibration peak, which represents the uncertainty in the 
measurements due to drifts, usually amounted to less than 0-25%. 
Spectra of K x-rays from natural lithium, *Be, 1B (95%) and 1'B 
(82°) are shown in figs. 3-6}. Kg peaks are seen clearly resolved in the 
spectra of lithium and beryllium x-rays. The boron spectra, taken with 
a counter filled to 5 atmospheres pressure of xenon to obtain an increased 
efficiency, also contain peaks due to boron » K x-rays and aluminium 
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M x-rays from the material of the source box. For elements of higher 
atomic number, the yield of x-rays and the efficiency of the counter 
were too low for worthwhile measurements to be made. Peak positions 
were determined from free hand drawings through the experimental 
points by two observers, the mean position being used. Only the K, 


Fig. 5 
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peaks were measured as there is evidence of unresolved higher transitions 
in the Kg peaks. Table 2 lists the energy values obtained. The mean 
value quoted is the mean of the individual values weighted by the number 
of counts in the peak in each case. The errors given are estimated 
ones based on the uncertainty in the peak position, and random drifts. 
They are somewhat larger than the r.m.s. errors derived from the internal 
consistency of the results. Values of the K., energies of individual isotopes 
are derived from the measurements by means of the relation 
H=a,H,+a,H, 

where £ is the energy observed with the mixture of isotopes, H, is the 
energy corresponding to the isotope having a percentage abundance @,, 
and similarly for #, and a,. In each case the correction due to the 
admixture of the second isotope is very small. 

In table 3 the derived K, energies for individual isotopes are given 
together with the values obtained by Stearns and Stearns (1956 a) and, 
the theoretical values (m,—272-8 m,) neglecting the short range pion 
nucleon interaction but including finite nuclear size and vacuum polariza- 
tion corrections. The agreement with Stearns’ values is excellent except 
in the case of ''B and confirms the existence of an energy shift due to 
the short range pion nucleon interaction. As can be seen from fig. 6 
there are no signs on the proportional counter spectra of any additional 
lines above the K, peak which could affect the measured energy in the 
case of 4B, The sample used was of 99-5°,, purity and showed no signs 
of a peak at 74 kev due to Mg L, rays, which was prominent when less 
pure samples were employed. 

The energy shift, which is due to short range interactions between the 
pion and the nucleus, has been calculated by Deser et al. (1954) in terms 
of the interaction of the pion with single nucleons summed over all the 
nucleons in the nucleous. 

Brueckner (1955), considered in addition the interaction of the pion 
with two nucleons, the process leading to star formation, and predicted 
energy shifts approximately 1-5 times as large. Karplus (1956), pointed 
out that the additional shift predicted by Brueckner is very sensitive to 
the momentum distribution of nucleons within the nucleus. He favoured 
a distribution giving a smaller shift than Brueckner’s. 

The percentage shift in energy given by Deser et al. is 

9 Ve 
oe = = : = 43 Za,+ es Fas} 
Ke 3d Ty |B 3 
where a, and a, are the pion nucleon s scattering lengths for isotopic spin 
states 1/2 and 3/2 respectively, 7, is the mesonic Bohr radius and p and 
are the reduced masses of the pion nucleus and pion single nucleon 
systems. 

The percentages energy shifts given by this equation, using the values 
a,=0:167--0-012 A, and a,=—0-105+0-010A, (Orear 1956), where A, is 
the 7-meson Compton wave-length, are given in table 4 together with the 
experimental values. The calculated shift depends very nearly on 
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Table 2. 


Element and 
isotopic 
constitution 


Lithium 
Enriched in 
lithium 6 


Lithium 
("Li 92:5%) 


(ola 7-647,) 


Beryllium 
*Be 100%) 


Boron 
(B 95%) 
(MB 5%) 


Boron 
fB Sl-6%,) 
296, 18497) 


Energy 
calibration 
used 


Measured values 
of w-mesonic 
K, energy (kev) 


24-13 kev} 


(2488n) 


24-13 kevt 
(138) 


41-31 kev 
(153Gd) 


2-19 kev 
(20H g) 


72:19 kev 
(Hg) 


4 determinations 


23-Be 
23-9, 


Measured Energies of 7-Mesonic K, Lines 


Mean 
observed 
K,, energy (kev) 


23°88-+0-1 


far) 
Or 
We 
~l 


| (-4. 


u 


64-65 +0°3 


+ A critical absorber (silver) removed the Ky radiation. 
t Measurement of the escape peak. 


Table 3. 
Measured K,, 
Isotope | energy (kev) 
present 
work 


(Stearns and Stearns 1956 a) 


Measured and Calculated Knergies of 7-Mesonic K, Lines 


Measured K,, 
energy (kev) 


Calculated K, 
energy (kev) 


(E.M. effects only) 


Mz=2712°8 mM, 


SLi 24-31 +0-1 
“le 23-85 +0:1 


’Be 42-15+0-15 
cad b3 65-3. +0-4 
75 64:5 +0°3 


23-77 +0-12 
42-09 +0-1 
65:2 +0-2 
63-5 +02 


24-53 
24-61 
43-95 
68-75 
68-82 
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a,+2a, and thus involves a small difference between two quantities. 
The resulting errors in the calculated shifts are considerabley. te 
Brueckner theory predicts an additional percentage shift Z?/985 . 4/3 
and the corresponding values, together with those of Karplus and. 
Halpern (1957), given in table 4 are subject to a similar uncertainty. 


Table 4. Measured and Calculated Energy Shifts 


Observed energy Se ans Calculated 


Calculated | Calculated 


Tsotope shits Stearns and ae Brueckner | Karplust 
present results Stearns (1956) et al. 

7o % % % 70 
6Li 0-9-L0-4 1140-5 | 2340-5 | 1-4+40-5 
"hi 3-140-4 3-4-£0°5 2-540:7 | 3-7-0-7 | 3-1-£0-7 
*Be 4-1£0:3, 4:24-0:25 3941-1 | 6241-1 | 4-941-1 
0B 5-0-40-6 5-1540°3 3241-6 | 6641-6 | 39416 
1B 6-3-£0-46 7-8+40°3 5541-7 | 8941-7 | 6941-7 


t Taken as 1.25 times the values of Deser et al. 


Table 5. Comparison of the Theories 


a,+2-23a3 a,+2-13a, 
Theory Prop. counter Stearns and Stearns 
results scintillation counter results 
Deser et al. —0:079 +0-004A,, —0-069 +0-002A, 
Karplus —0-063 +0-003A, —0-055 +0-002A,, 
Brueckner —0-033 +0-005A, —0-023 +0-002A, 
Orear’s value —0:067 +0-0234A, —0-057 +0-022A, 


The uncertainties in the calculated values of the energy shifts makes a 
clearcut discrimination between the theories impossible. 

To illustrate this the experimental shifts from all the isotopes were 
added and a value for a,+2-23a, was deduced according to the three 
theories. These values are given in table 5 together with the value of 
the same quantity derived from Orear’s values of a, andaj. Stearns’ 
data, which includes measurements up to Z=9, gives values of a,-++2-13a, 
also included in table 5. 

Thus both the theories of Deser et al. and Karplus are equally in 
accord with the experimental facts, and Brueckner’s theory is not com- 
pletely excluded although it gives the worst fit to the data. 

Stearns used their data to deduce values of a, and ag assuming that the 
theory of Deser et al. applied. A least squares fit gave values of 
a,=0-145A, and a;=~0-10A, which agree well with Orear’s values. It 


+ Account was taken of the fact that the errors on a, and a, given by 
Orear are not independent. We are indebted to Professor H. L. Anderson for 
pointing out the possible importance of this and to Dr. J. Howlett for 
instruction in the method of combining such errors. 
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does not however follow, as is evident from table 5, that the theory of 
Deser et al. gives the best fit to the data. When the errors in the least 
squares values of a, and as (+0-02 and --0-014 respectively according 
to Orear 1956) are taken into account, the values of the energy shifts 
deduced are seen to be subject to a very large error. For example the 
calculated shift for oxygen is 10-4+~ 6%. 

The present experimental data is probably accurate enough to discri- 
minate betwen the various theories when more precise values of a, and 
a, become available. 


§ 5. ENERGIES OF 7-MeEsonic L X-RAYS 

The effect of short range forces on the energies of L transitions from 
light elements is very much less than in the case of the K series. Stearns 
and Stearns (1956 a) have estimated that an upper limit to the energy 
shift is given by —éH#/E,,=2:5x10->xZ*%. The important phase 
shift (i.e. for the 2p level) is 5,, which is positive so that the energies of 
the L transition would be raised in contrast to the behaviour for the 
K transitions. Stearns and Stearns (1956 a) found no evidence for a 
shift of this magnitude for Ca(Z=20), and it is therefore safe to assume 
that the actual shift is much less than the upper limit. For Z up to 
about 10 the L ray energies can therefore be calculated directly from 
electromagnetic effects to better than 0-1%. 


Fig. 7 
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L series x-rays from carbon. 


Energy measurements were made of the L, (3— 2) lines of elements 
from carbon to magnesium, primarily as a check on the experimental 
technique which was similar to that for the K series measurements. 
The yield of L rays is however higher than that of the K rays and shorter 
runs were possible—four hours being about the maximum period. 

L ray spectra from C, N, Na and Mg are shown in figs. 7-10. The 
Lg(4— 2) and in some cases the L, (5 — 2) groups are resolved. In the 
cases of sodium and magnesium the M,, (4-3) line is also clearly seen, 


Fig. 8 
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L series x-rays from nitrogen (C,N,H,). 
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L and M series x-rays from sodium. 


Fig. 10 
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L and M series x-rays from magnesium, 
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Table 6 contains the measured L, energies together with the values 
calculated for m,—272-8 including the small correction due to vacuum 
polarization (Stearns and Stearns 1956 a). The finite nuclear size 
correction is negligible. From the mean value of the measured energy 
an effective meson mass is deduced from the measurements for each 
element. The values are given in the last column of Table 6. 


Table 6. Measurements of L,energies 


Measured Calculated | Etecttve 
energy values Weighted 7-meson 
Element | Compound | ~ y. . | mean energy ae mass 
of L, line Mx =272°8 
(kev) — (kev) (units of 
Me) 
18-1, 
C 18-1, 1S. 18-38 269-4 
3-17 
BN 24-9, 
NH,F 24-6, 
N NH,F 24-8, 24-8, 25-07 270-4 
BN 24-9, 
Lif 41-6, 
Lif 41-9, 
NaF 41-9, 
NH,F 41-6, 41-7, 41-60 273-9 
F NH,F 41-6, 
Lif 41-7 
LiF 41-7,+ 
LiF 41-7,7 
61-6, 
Na 62:1, 62-0, 62-29 271-9 
62-2, 
74-2 . 
Mg 74-0, 74-1, 74-18 272-7 


+ Measurements of the ‘ escape’ peak. 


The mean of all these values weighted by the total number of counts in 
the peaks divided by the instrumental width squared is m,—=272-5+-0-8 m, 
which is to be compared with the accepted value 272-8-+-0-3 m, (Barkas, 
et al. 1956) and the limits 272-2 m,<m,,<273-6 m, set by critical absorber 
measurements on z-mesonic M ray lines by Stearns, Stearns, DeBenedetti 
and Leipuner (1954). 

The value of m, deduced from the L ray measurements therefore 
agrees with the accepted value within the experimental error, It 
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establishes that there is no large systematic error in our technique of 
measuring energy. In addition there is no indication within the accuracy 
(--0-6%) of any shift in the energies of the L rays due to short range 
interactions up to Z=12. 


§ 6. MEASUREMENT OF THE NATURAL BROADENING OF THE a-MESONIC 
K, Live From BERYLLIUM 


The strong interaction between 7-mesons and nucleons which leads to 
capture from the Is or higher levels is also expected to broaden the 
energy levels, especially the 1s level. Brueckner (1955), has calculated 
that the 1s level should have a full width at half maximum 
P=2~ 4/3 x Z?/2150. Hy,, which gives [=0-9 kev for beryllium. The 
instrumental width in our proportional counter was equivalent to a full 
width at half maximum of 3-4 kev at this energy. However the effective 
instrumental width can be reduced by making use of the escape peak. 
Here the instrumental width is only 1-6 kev at the energy (12-5 kev) 
corresponding to the escape peak from the beryllium K, line. Natural 
broadening of the line will be equally great for the main and escape 
peaks since there is a constant difference of energy (equal to that of the 
electronic K, x-ray of xenon) between the two peaks. The instrumental 
width is therefore small enough to enable a natural broadening of the 
predicted amount to be detected. 

It was desirable to measure the instrumental width in conditions 
approximating as closely as possible to those in which the beryllium 
K, line was measured. Accordingly a a-mesonic L line (fluorine 
L,=41-6 kev) of very nearly the same energy as the beryllium K line 
(42-15 kev) was chosen for this purpose. The expected natural broadening 
of this L series line is negligible and the observed width is essentially 
instrumental. 

Beryllium and lithium fluoride sources of very nearly the same mass/cm?2 
and area were used, and the proportional counter single rates were made 
as nearly as possible the same for the two runs. Two hour runs with 
beryllium alternated with one hour runs with LiF, with calibrations 
between each run, and the individual spectra were added at the end of a 
series. Hence drifts should affect the two spectra equally. In fig. 11 
two spectra taken in this manner are shown. Figure 12 shows similar 
spectra in which only the escape peaks were measured. 

The widths of these lines were measured by fitting Gaussians to the 
spectra using an ogive plot, on which the fractional area under a Gaussian 
up to any abscissa when plotted against the abscissa is a straight line. 
The fit is sensitive to the background assumed and _ several choices 
indicated in figs. 11 and 12 were made. 

In table 7 the measured values of the standard deviations o obtained 
from the ogive plots are given. In addition the standard deviations of 
the calibrations (41-3 kev and 11-6 kev for the main and escape peak 
measurements respectively), taken between the individual runs are 
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given. They represent the broadening due to random drifts in gain and 
apply equally to the Be and F measurements. The drift is in all cases 
negligibly small. The natural broadening of the beryllium line onat 
is then given by onat=(op.2—c . o,2)1/2 where the factor c is 1-01 for 
the main peak measurement and c=1-04 for the escape peak measure- 
ment and takes account of the slight difference in energy of the two 


Fig. 11 
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Comparison of the widths of the main peaks from B,, Ka and F, L, radiations. 


i 


lines. The experimental values are given in the last column of table 7 
with errors assigned by combining the largest value of the beryllium 
width with the smallest value of the fluorine width and wice versa. 
The weighted mean of these values is o=0-52-0-06 kev ort 
T=2-36c=1-2--0-2 kev which is close to Brueckner’s value [’=0-9 kev. 


+ Ina previous letter West and Bradley (1956), J” was wrongly equated to o, 
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The experimental width corresponds to a mean life of 6x 10-1 sec for a 
7-meson in the 1s state of beryllium. 

The Brueckner theory, therefore, which somewhat over-estimates the 
energy shift, is nevertheless in agreement with experiment as far as the 
width of the line is concerned. A possible explanation pointed out by 
Eden (private communication) is that the width as calculated by Brueckner 
is much less sensitive to the momentum distribution of nucleons than is 


Fig. 12 
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Comparison of the widths of the escape peaks from Be, K, and F, L,, radiations. 


his energy shift. No level broadening is expected from the pion, single 
nucleon interaction responsible for the level shifts calculated by Deser 
etal. (1954) since the charge exchange process which could give a broadening 
is energetically impossible in °Be. 

The fact that nuclear absorption of 7-mesons does occur and moreover 
that the resulting broadening calculated by Brueckner is in agreement 
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Fig. 13 
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Ogive plots of the escape Peaks from fluorine L, and beryllium K, Radiations. 


Table 7. Natural broadening of the Be K, line 


. o 7 2\1/2 
See Calibrations | °Be oF | Onat= (eRe eo") 
details (kev) (kev) (kev) (kev) 
Main peak 
measurement 0-2 1-8, 1-44 
Be, two 2 hr runs 1-7, +1:76 1-55 >1-46 1-0 -+0-2 
F, three 1 hr runs 1-63 1-44 — (0-4 
Escape peak 0:8, 0-68 0-69 0:47 -+-0-10 
measurement 0-8, >0-85 0-7, ‘ 
Be, four 2 hr runs 0-08 0-8, 
F, three 1 hr runs 
Escape peak 0:86 OO aGn. a 0:49 -+-0-07 
measurement ps 0:8, >0:84 0-65 0-66 
Be, four 2 hr runs 0-08 0-83 


F, four 1 hr runs 
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with experiment makes it unlikely that the theory of Deser et al. can 
adequately describe the energy shifts of the K lines. As discussed in 
§ 4, the energy shifts themselves cannot decide between the theories at 
present but the existence of a level broadening is evidence that a theory 
of the Brueckner or Karplus type which includes the interaction leading to 
absorption of the z-meson is to be preferred. 


§ 7. RELATIVE INTENSITY MEASUREMENTS 


The proportional counter is not a suitable instrument for the absolute 
measurement of intensities of the x-rays owing to the rapid variation of 
efficiency with energy and the very complicated geometry. Relative 
efficiencies for two lines of slightly different energy can however be 
calculated with fair precision from the photoelectric absorption co- 
efficient of xenon as a function of energy. Only an approximate value 
of the mean track length of x-rays passing through the counter is required. 
It was derived from a measurement of the efficiency at 60 kev made 
using an 744Am source of known intensity placed at the same distance 
from the counter as the mid point of the target. The mean track length 
determined in these conditions was 0-7 (counter diameter) for X-rays 
crossing the counter. 


Table 8. Relative Intensities 


Element Intensity ratio Intensity ratio 
i ° Ls aE. 


ba 


1. 0:35 
We 0-30 
mean 0:32-+ 


Beryllium 


mean 0-20; -+0-02 
Carbon ll. 0:75 : 0-20: 
Oe 0:89 : 0:09 : 0: 
mean 0-82-+0-08 : 0:14+0-05 : 0-03 


Nitrogen i OTT OT 
Op 0°83 20-13% 
mean 0-80-40-05 : 0-15 -+0-02 : 0-05 4 


Absorption in the target material, the anticoincidence scintillator and 
the aluminium wall of the counter was approximated by assuming first 
that the actual thickness of the scintillator and aluminium and half the 
thickness of the target was effective. A second calculation was made 
using twice these values to allow for obliquity. 
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Table 8 contains relative intensities of K, (2—> 1) and Kg (8 > 1 and 
4—> 1 ete.),-and of L,, La, L, lines deduced in this fashion where 1 and 
2 refer to the two assumptions made in correcting for absorption. 

The Kg: K,+Kg ratio in beryllium is similar to the values observed 
in nitrogen and oxygen (0-20-+0-05 and 0:21-+0-06, Camae et al. (1955)), 
but the ratio for lithium is definitely greater than for the heavier elements. 
This may be connected with the fact, observed by Stearns and Stearns 
(1956 b) that the yield of lithium K series x-rays is smaller than that 
from beryllium, which they ascribe to Auger competition. 

The relative yields of the L series lines agree with the scintillation 
measurements of Camac, Halbert and Platt (1955). 

A comparison was also made of the total K series yields of ®Li and 7Li. 
The triple coincidence rates were monitored for these runs which were 
interlaced with each other and the number of mesons stopping was 
assumed to be proportional to the number of electrons in each target, 
and to the number of mesons passing through counters 1, 2 and 3. A 
ratio 

= StH from “Li per stopped meson —~0-93-4.0-05 
K x-rays trom ‘Li per stopped meson 


was obtained. 
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ABSTRACT 


It is shown that the rotation axis of a boundary containing dislocations 
of three systems can lie anywhere on the surface of a cone; to each 
direction of the rotation axis there corresponds a single boundary plane. 
The theory cannot explain the experimental evidence on misorientations 
in zinc unless slip occurs in non-basal directions, or on planes not hitherto 
reported, at high temperatures. The experimental evidence for face- 
centred cubic metals is insufficiently precise to form a test of the theory. 


§ 1. [INTRODUCTION 


In a previous paper (Ball and Hirsch 1955) it has been shown that Frank’s 
equation for the dislocation content of a plane, low angle boundary 
(Frank 1950) can be used to determine the types of stress free boundary 
that can be formed with dislocations of one or two systems subject to 
these dislocations lying in their usual slip planes. The rotation axes 
of most of the boundaries observed experimentally approximate fairly 
closely to those of boundaries containing only two sets of dislocations, 
but often the discrepancy is greater than can be accounted for by 
experimental error. In this paper the method of analysis is extended to 
boundaries containing three sets of dislocations, and the theory is applied 
to the experimental observations. Some new observations are reported. 


§ 2. METHOD of ANALYSIS 
Frank has shown that the resultant Burgers vector d of the dislocations 
cut by a vector r lying in a stress free boundary is given by 


d=(r X1)2 sin $a 


where I is a unit vector parallel to the axis of rotation and « is the angle 
of rotation. It follows that dis normal to r and I for all r. This relation 
can be used to determine the types of boundary that can be formed with 


+ Communicated by the Author. 
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any combination of dislocation systems. The most important cases 
involving three sets of dislocations will now be discussed. 

Let the normals to the three slip planes be p;, Ps, Ps, the Burgers 
vectors b,, by, bs, the dislocations parallel to vectors rj, fp, rs and let 
there be 1, %», Nz dislocations of each type per unit length taken normal 
to that type. 


2.1. Two Slip Systems have the same Slip Plane 


Let the identical planes be p;, p». If we take r parallel to r,, it will 
cut dislocations of the third set only, so r,; must be normal to bs. By 
taking r parallel to rz, we see that | must also be normal to the vector 
(n,b,+7.b,). As this vector can lie anywhere in the plane containing 
b, and b,, | can lie anywhere in the plane normal to b,. 

Corresponding to each direction of the rotation axis there is only one 
possible boundary plane; for r,; must be normal to b, and lie in p,, 
and r,; must be normal to (n,b,+7,b,) and lie in p;. The boundary 
plane is the plane containing r, and rg. 

If the dislocations are all parallel, they must be parallel to the inter- 
section of the planes p,, p;. | will also be parallel to this direction. The 
boundary plane must be normal to the vector (n,b,-+-n,b,+nsb,). If 
b,, b,, b, are coplanar, the boundary can only form on one plane. Other- 
wise, any boundary plane is possible that is parallel to rg. 


2.2. Two Slip Systems have the same Burgers Vector 
Let the identical Burgers vector be b,, by. Then I must be normal to 
the plane containing b, and bg, as d will always lie in this plane. 
r; is determined as it must be normal to b, and lie in plane ps. The 


boundary plane must contain rz, but otherwise it is arbitrary as n, and n, 
are independent. 


2.3. The Three Slip Systems are Unrelated 
It will be convenient to define unit vectors d,, dy and ds by the relations 


with similar equations for d, and d, obtained by rotating subscripts. 

If now r is taken parallel to r,, d will be parallel to d,, as it must lie 
in the plane containing b, and bs and be normal to I. 

Hence r, can be determined, as it must be normal to d, and lie in the 
plane p;. Similarly, we can determine r, and rz. If ry, r, and rs so 
obtained are coplanar, the boundary will be satisfactory. 

In practice, r,, r, and r, are obtained in terms of the crystallographic 
indices of the rotation axis (7, K, L). The indices of the plane containing 
r, and r, are next found, again in terms of H, K, L, and the condition 
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that r; should lie in this plane gives an equation relating H, K and L. 
Solving this equation gives the locus of the rotation axis. The method 
fails for certain directions of the rotation axis, which must be checked. 

It is not immediately obvious that the dislocation ratios required by 
the above construction will be mutually compatible. That they will be 
can be shown thus. Read (1953) has shown that any boundary can 
be composed of any three sets of dislocations with non-coplanar Burgers 
vectors, the direction of the dislocations being arbitrary. A boundary 
obtained by the above construction will require its dislocations to lie 
in their usual slip planes, for, taking r parallel to the intersection of the 
boundary plane and plane p,, d will be parallel to d,, and so r must be 
parallel to r, as it cuts no dislocations of set 1, i.e. r, lies in the plane py. 
Similarly r, and r; must lie in the planes py, pz respectively. 


§ 3. RESULTS 
3.1. Hexagonal and Body-Centred Cubic Systems 


The only slip directions for hexagonal metals reported in the literature 
are of the form ¢<1120) lying in the basal plane. The most common slip 
plane is (0001) but slip has also been reported on prismatic and pyramidal 
planes, particularly at high temperatures. The only boundaries that 
can be formed with dislocations of these systems are (a) twist boundaries 
with rotation axis [0001] or (6) tilt boundaries in which all the dislocations 
are parallel. The rotation axis for such a boundary is parallel to the 
intersection of the slip planes, while the boundary plane contains this 
direction and is normal to the basal plane. 

In body-centred cubic metals slip takes place in the direction [111] 
on a variety of planes. The number of possible combinations of slip 
systems is very large, and they will not be considered in detail. If the 
slip planes are suitably chosen, practically any boundary can be formed 
by slip alone. 


3.2. Face-Centred Cubic System 


In face centred cubic metals slip takes place on {111} planes in (110) 
directions. We will now enumerate the distinct combinations of slip 
systems, using Thompson’s notation (Thompson 1953). 

There are five distinct combinations of vectors 


1. Three vectors equally inclined to one another. 
2. Two vectors at right angles and a third. 

3. Two vectors identical and one at 60° to them. 
4, Two vectors identical and one at 90° to them. 


. Three vectors in a plane. 


Or 
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Combinations of vector and plane are 


a b c d 
Case | A B y y 
LAE 5) B 
aD B B 
Case 2 AVS 3) y 
BC 0 o 
AD p y Y 
Case 3 A B y Y 
AB 8 rs) 
AD B 
Case 4 A B y 
AB 8 
CD ¢. 
Case 5 AB 5 ) r) y 
BeC 5) ) oa io 
GLA 8 B B B 


Most of these combinations are of the special types considered in § 2.2 
and § 2.3. 
Possible rotation axes are 


a In (110) Cases 1b, 2c 

b [111] Cases 3a, 3b, 5c-d 

c [100] Case 4 

d In (111) Case 5a 

d [110] Cases 1b, 2c, 3b, 5b (pure tilt boundaries). 


Cases la, 2a, and 2b are of the type considered in § 2.4. The loci 
of the rotation axis for these cases are shown in figs. 1, 2 and 3, respectively. 
The equations to be satisfied by H, K, L, are 


la. H?4K3418—2H?L—2K*L—21°H—HKL=0 
2a, 3L§+H%(K+L)—K*X{L+H)+2L°(2H+K)+3HKL=0 
2b. K8—218+ H*(K+L)+2K%(H+2L)+L°H+K)+3HKL=0, 


(The orientation of Thompson’s tetrahedron adopted was «=(111); 
B=(111), y=(111), 5=(111)). 
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Locus of rotation axis for case 2a of the face-centred cubic system. 


To each direction of the rotation axis, there corresponds a single 
boundary plane. This can be determined graphically, as in fig. 4, or by 
substituting values of H, K, L into the indices of the plane containing 
r, and rp. 


Locus of rotation axis for case 2b of the face-centred cubic system. 


Fig. 4 


Graphical determination of boundary plane. 


Fig..5 


Variation of densities of dislocations along the locus of the rotation axes 
from [111] to [111] for a boundary of type la. 


Surface Distributions of Dislocations in Metals : II 983 


The numbers of dislocations of each type in the boundary can be 
determined analytically (Read 1953) or graphically. Figure 5 shows 
the way in which the densities of the three systems of dislocations in a 
boundary of type la vary as the rotation axis travels the locus from [111] 
to [111]. It will be seen that boundaries containing only one or two 
sets of dislocations are included as special cases. 


Fig. 6 


Observed rotation axes in polycrystalline aluminium deformed at temperatures 
in the range 0-300°C. 


§ 4. Discussion 


4.1. Hexagonal and Body-Centred Cubic Systems 


There is very little experimental evidence on the rotation axes of 
low angle boundaries formed by deformation. The most extensive 
observations are those of Cahn et a/. (1954) on zinc. In polycrystalline 
zinc made to creep at 250°C most of the rotation axes observed were 
within 10° of the basal plane and it has been suggested (Ball and Hirsch 
1955) that this may be because the dislocations are heavily jogged. 
However, in specimens made to creep at 350°o the rotation axes observed 
were distributed fairly uniformly throughout the unit triangle. These 
boundaries could not have been formed by slip alone on any combination 
of known slip systems. They could have been formed with the help 
of climb, but the amount of climb required is very large and it seems 
more reasonable to suppose that at these elevated temperatures slip 
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is occurring in some direction that does not lie in the basal plane, or 
on some hitherto unreported planes, 

All observed rotation axes in body-centred cubic metals can be ex- 
plained in terms of boundaries containing only one or two sets of 
dislocations. 


4.2. Face-Centred Cubic System 


The rotation axes shown in fig. 6 have been observed in polycrystalline 
aluminium deformed at temperatures in the range 0-300°c. Details of 
these experiments will be given later. 

It will be seen that most of the rotation axes correspond fairly closely 
with those of boundaries that could be formed with dislocations of only 
two systems ([100], [111], in (110), in (111)) but that in some cases the 
discrepancy is greater than can be accounted for by experimental error. 
It is thought that most of these boundaries contain a small admixture 
of dislocations of a third system. Boundaries of types la, 2a, 2b (§ 3.2) 
give such complete coverage of the unit triangle that, with the precision 
of the present experimental results, any rotation axis could be explained 
in terms of one at least of these types of boundary. However, as to 
each direction of the rotation axis there corresponds only a single boundary 
plane, a test of the theory could be made if observations were possible 
of both rotation axes and boundary planes. This is being attempted. 
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ABSTRACT 


Using a generalization of the model first introduced by Mott for 
dealing with impurities in monovalent metals, the interaction between 
an impurity and a vacancy has been studied. Detailed computations 
have been carried out for the case of copper, with a divalent impurity 
and a vacancy at a separation of 5 atomic units. The interaction is 
attractive, the associated energy being calculated as 0-08(4) ev when 
exchange and correlation energies are neglected. Inclusion of these 
effects reduces the interaction energy to 0:03(3) ev. A previous calcula- 
tion by Lazarus leads to a value of 0-045 ev for the case considered here, 
while a modification of the Lazarus expression proposed by the present 
writers gives 0-035ev. Thus the present work tends to substantiate 
further the Lazarus theory of solute diffusion in metals. 


§ 1. LyTRODUCTION 


Mucu experimental information has accumulated over the past few years 
concerning the diffusion of impurity atoms in metals. Particularly 
interesting from the point of view of the present work is the finding 
that the activation energy for diffusion of cadmium, indium, tin and 
antimony in silver decreases as the valence difference between solute and 
solvent atoms increases (Sonder ef al. 1954, Tomizuka and Slifkin 1954), 
Later work by Tomizuka and Sonder (1956) on self-diffusion in silver 
also contains relevant results and should therefore be mentioned here. 

An important step towards an understanding of the salient features 
of these experimental results was taken by Lazarus (1954; see also 
Blatt 1955, Alfred and March 1956) who proposed an explanation based 
on the screening of a solute atom in a metal. Fundamental to the 
Lazarus theory is the evaluation of the interaction energy between the 
screened solute atom and a vacancy, and a simple electrostatic argument 
was used to compute this. However, as Lazarus pointed out elsewhere 
(Lazarus 1955), such a procedure appears to neglect any changes in 
kinetic energy and also to neglect interaction between the screens of the 
solute ion and the vacancy. (see, however, § 3), as well as refinements 
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such as exchange and correlation effects. In spite of these and other 
approximations which it was necessary to make, Lazarus’ theory has 
had a good deal of success. In addition to correctly explaining the 
variation of activation energy for impurity diffusion with valence difference 
referred to earlier, the theory also predicts a correlation between the 
frequency factor D, and the activation energy H in the diffusion equation 


D=D, exp (—H/RT) sie a 


where D is the measured diffusion coefficient, R is the gas constant and 
T is the absolute temperature. Such a correlation is indeed evident from 
the experimental results and for diffusion in silver reasonable quantitative 
agreement between theory and experiment is found. Even for diffusion 
in nickel, to which the present theory can hardly be quantitatively 
applicable, such a correlation is again found and the agreement between 
theory and experiment is better than order of magnitude (Swalin 1956 ; 
see also March 1957, where a short review of this work can be found). 

In view of these successes, it seemed a matter of some importance to 
attempt to put the theory on a sounder basis by making a more careful 
calculation of the interaction between a solute atom and a vacancy in 
ametal. For convenience later, we note first of all that for the interaction 
energy 4H between an impurity having Z+-1 electrons outside a closed 
shell and a vacancy in a monovalent metal, with Fermi energy £,,, 


Lazarus takes 
Z 25/2 
AH=— —exp (—q’) ; Ga — hy eee) 

re 7 
7) being the impurity—vacancy separation. This is obtained by calculating 
the screened field of the impurity from the linearized Thomas—Fermi 
equation (see, for example, the review by Friedel 1954). Blatt (1955) 
attempted to improve the treatment by regarding the screening constant ¢ 
as a parameter to be determined by the requirement that the total 
screening charge should be correctly given, following an earlier proposal 
made by Friedel (1954). Unfortunately, much of the agreement between 
theory and experiment then disappeared, the variation of the activation 
energy H with Z being much too small. However, the present writers 
(Alfred and March 1956) subsequently showed that, if the exact solutions 
of the Thomas—Fermi equation for an impurity are used (Alfred and 
March 1955; see also Fujiwara 1955, who independently made some 
similar calculations) instead of Blatt’s fields, then the interaction energy 
may be written 
Z 

— me exp (—q1o) eee TS me es) 
where « varies with Z, but is of the order of unity ; for silver with Z=1 
for example it was found that ~=0-75. Using these results, the agreement 


+ In (2) and all subsequent equations, atomic units are used, 
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between theory and experiment was largely restored although, as the 
present writers pointed out, the formula of Lazarus is probably preferable 
on the grounds of internal consistency. 

From the above discussion, it is clear that ideally we should like to 
know the interaction energy as a function of impurity—vacancy separation 
around the nearest-neighbour distance in a metal, but the calculations 
turn out to be formidable and we have been obliged to restrict ourselves 
to the more modest programme of obtaining the interaction energy for 
one separation, which we have taken conveniently as 5 atomic units. 
This is close to the nearest-neighbour distance in copper (actually 
4-8 atomic units) which is the metal we consider in our computations. 

In addition to the work referred to above, other attempts have been 
made previously to calculate the interaction energy between imperfections 
in metals. For example, we might mention the work of Bartlett and 
Dienes (1953) in which the association energy of two vacancies in copper 
was estimated to lie between 0-2 and 0-6 ev. More recently Seeger and 
Bross (1956) have made a much more detailed investigation of this 
problem and find a value around 0:3 ev. 


§ 2. MopEL FoR PRESENT CALCULATION 


We shall now outline briefly the model adopted for the present 
calculation. We shall deal with the case when the solvent is a monovalent 
metal, and shall describe the conduction electrons in the unperturbed 
metal by the free electron model, neglecting exchange and correlation. 
However, we shall eventually consider the effect of these energy terms 
by a less accurate, but probably adequate, method. Following Mott 
(1936) we represent an impurity with (Z+1) electrons outside a closed 
shell as a point singularity carrying charge Z. As Mott showed, and as is 
essentially embodied in the statement of eqn. (2), a self-consistent field 
calculation then leads to the result that the conduction electrons of the 
metal adjust themselves to screen out the field of the charge Z in a 
distance of the order of 1A. For the vacancy we shall adopt a model 
which appears to have been first suggested by Dexter (1952) ; namely 
we shall treat the vacancy as an impurity of zero valency. Thus for the 
monovalent metal we are considering here, the vacancy is represented 
as a point singularity carrying a single negative charge. 

To obtain a first approximation to the solution of the present problem, 
we need the exact numerical solutions for a single impurity and a single 
vacancy in an infinite metal. As we have already remarked, the case 
of an isolated impurity has been treated by means of the Thomas—Fermi 
approximation, first used in this connection by Mott, and a number of 
exact numerical solutions are available for copper and silver (Alfred 
and March 1955, Fujiwara 1955). Unfortunately the exact numerical 
solutions referred to above are all for the case of positive Z. For negative 
Z, the present writers (Alfred and March 1956) have given a very simple 
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analytical argument by means of which solutions of fair accuracy can 
be obtained in this case and this can be made the basis for the calculations 
of the exact numerical solution (see § 4). Thus the potentials round an 
isolated impurity and an isolated vacancy, along with the corresponding 
electronic charge distributions, were known at the outset. 

The major problem, therefore, consists in the redetermination of the 
distribution of conduction electrons and the modified self-consistent field 
when the impurity (we shall deal specifically with the case of a divalent 
impurity, say zinc, for the purpose of the numerical calculations, although 
the method adopted is obviously generally applicable for any valency) 
and the vacancy are separated by a distance of 5 atomic units. This is 
a problem which from a mathematical point of view at least is not entirely 
dissimilar to that of determining the distribution of electrons in a hetero- 
nuclear diatomic molecule, and we were encouraged to attempt the 
numerical solution of the present problem as a treatment of the iodine 
chloride molecule had previously been carried through successfully using 
the relaxation method (Glazer and Reiss 1953). We have, therefore, 
made use of this technique in our work, even though an additional 
complication, that of a variable boundary round the vacancy (see § 5), 
was present in our problem. 


§ 3. ImpuRITY—VaAcANcYy INTERACTION ENERGY 
IN Frrst-ORDER TREATMENT 


In view of the success of the simple formula (2), it is of interest at this 
stage to consider a more detailed calculation of the interaction energy 
between an impurity and a vacancy using the model outlined above, but 
employing the Thomas—Fermi method only in the linearization or first- 
order approximation (Mott 1936; see also Friedel 1954). In this case, 
heavy numerical computations are avoided and some insight into the 
main features of the problem can be obtained. On the other hand, there 
appear to be a number of objections to the first-order method, and we feel 
therefore that the quantitative predictions must be viewed with reserve. 
Dr. J. Friedel first brought most of the essential results obtained in this 
section to our notice, although we believe the detailed argument given 
below to be new. 

We take as the basic equation of the first-order treatment 

VAV = —4n(n—n = 0 Ve e e 
where JV, is the perturbing potential energy, n(r) is the electron density 
in the perturbed metal and ny is the unperturbed density. We now seek 
a variation principle which will yield the relation (4) between density 
and potential. The form of this is suggested by the following argument. 
In the full Thomas—Fermi approximation the change in kinetic energy 
of the electrons is given by 

3 


0 (30228 | (7 2!3— 7,912) dr. 4 ; : , : : (5) 


— 
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Assuming 4=n—nz, is small, (5) may be written 
2 
Hy, | A det = | A* de Rr Be hei) 


with neglect of higher terms. Adding the potential energy change to (6) 
and requiring that the resulting expression is stationary with respect to 
variations in 4 then yields exactly the relation (4) between density and 
potential. Thus, we adopt (6) as the kinetic energy change consistent 
with the linearized scheme given by egqn. (4). 

It is now a relatively easy matter to calculate the interaction energy 
between two impurities in a Fermi gas. Suppose, to be general, the 
impurities are represented by point charges Z and Z’ in this model, 
the separation being 75. If V, is the perturbing potential due to Z alone 
and V,’ that due to Z’, then from (4) it follows that the potential in the 
case of two impurities at a finite separation is just the superposition 
potential V,+V,’. 

We must now compute the changes in potential and kinetic energies 
when we bring the impurities together, and we may write down the 
following contributions : 


(i) The interaction energy between the charge Z of one impurity and 
the perturbation V ,’ due to the other. 


(ii) The interaction energy between the charge cloud (—q?/47)V, 
round the impurity Z and the perturbation V,’ due to the other. 


(iii) The change in kinetic energy. 
These three terms are evidently given by 
ZL’ exp (—4"y) 
"9 


(i) 


2 
(ii) A a [V.V," dr 


oe 2a gq? ry , ' 
(iii) =| dal lo Aare {v2 i [v,rar}] 


qg : 
= 1 I{v,v, dr. 


Thus the contribution (ii) is just cancelled by the change in kinetic 

energy (iii), and we are left with the formula of Lazarus 

ZZ' exp (—qr 
ie LON ma (ay 
"o 

Unfortunately, however, if an attempt is made to include exchange 
along the lines of the argument given above, it turns out that the 
formula (8) is drastically changed. Nevertheless, as we shall see later, 
the simple arguments given here lead to results which have a number 
of features in common with those we find from our detailed computations. 


3X2 
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§ 4. PoTENTIAL AND DISTRIBUTION OF ELECTRONS RounpD A VACANCY 


We shall now attempt to use the full Thomas-Fermi method, avoiding 
altogether the first-order approximation used above, and as a necessary 
preliminary to the main problem we shall discuss in this section the 
results we have obtained for the case of an isolated vacancy. We wish 
to solve the equation 

7/2 

Eipiae = — (i= V2) sea) 
for the self-consistent potential energy V,. We note first of all that 
near the negative point charge, V, is large and positive and inside a 
sphere of a certain radius r,, will exceed the Fermi energy #,,. It then 
follows that, in the Thomas—Fermi approximation, no electrons are 
allowed inside this sphere, the region being one of negative kinetic energy. 
The radius r, is determined by the condition 


Vitti * ape 2 a ee 
and inside the sphere the appropriate equation to be solved is 
97/2 
Oita. ps2 
Vel oe TE BIAS 6 Ve eye ee 


The boundary conditions subject to which we must solve eqns. (9) and 
(11) are evidently 


V,790 as r>o fh OE ae) 


dV : 
V,, and a°5 continuous at r=r, 


Equation (11) can be solved immediately in this case and we have 


i) 95/2 
iui apa mye FERS Rat Sy bao (PS) 


where £ is a constant to be determined. On the other hand, eqn. (9) 
has to be solved numerically. We found a step-by-step numerical method 
convenient, trial and error being necessary to determine the correct 
value of 7, but fortunately the approximate solution considered previously 
(Alfred and March 1956) is a very useful starting point. We show the 
exact solution thus obtained in curve 1 of fig. 1 together with the 
approximate solution given earlier by the writers (curve 2) and the 
result obtained by a naive application of the first-order approximation 
(V,=1/r exp (—qr): see curve 3). Also shown in curve 4 of fig. 1 is 
the exact Thomas—Fermi perturbing potential for Z=1, while in fig. 2 


the electron densities round an isolated impurity and an isolated vacancy 
are shown. 


§ 5. SeLF-ConsistentT Fratp Rounp Impurrry—Vacancy System 


We must this time solve the basic equations for the potential energy V , 
assuming symmetry only about the axis joining the impurity and the 
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vacancy. As in the isolated vacancy case, we shall have a surface S 
surrounding the vacancy inside which there are no electrons, although 
we must now expect that it will be somewhat distorted from a sphere. 
Inside this surface we must solve eqn. (11) whereas outside we must 


deal with eqn. (9). The boundary conditions which the potential must 
now satisfy are 


V,.7>—— as 71> 0 


1 eer ae oe PE 
V ras as To > 0) ( ) 


Vi720 as 11> fe > Oj 


r|Vo| 


0.2 


0.1 


0 | Zz 3 4 
Perturbing potentials 


. rVp for vacancy—exact numerical solution. 

. *Vp for vacancy—approximate solution. 

2 fV p=exp (—9r). é 

r |Vp| for impurity (Z=1)—exact numerical solution. 


Rm Whe 
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Fig. 2 


n(r) 


0-04 


0.03 


0.02 


0.01 


0 | 2 3 


Electron densities round isolated imperfections. 
1. Impurity (Z=1). 
2. Vacancy. 


Fig. 3 


Notation for impurity—vacancy case. 
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together with the conditions that V,, and its first derivatives are 
continuous across S. The notation used is shown in fig. 3. 

As we have remarked earlier, we now assume for the purposes of the 
computations that the solvent metal is copper, that the impurity is 
divalent, and that the separation is 5 atomic units, or the approximate 
distance between nearest-neighbours in the lattice. 

We begin the solution for the perturbing potential V » corresponding 
to the impurity—vacancy case by removing the singularities at the 
impurity and the vacancy. This can be done by substituting 

pee eV mo ee oe (18) 
where V,; is the perturbing potential energy due to an isolated impurity, 
V,, that due to an isolated vacancy, and 4 is the change from the super- 


position potential energy. We may then write in the region outside the 


surface S 
27/2 
—[V?V,,+ V2V ,,— V3]= = (Em —V yi —V py +8)2?2—E p22] (16) 


where V?V,, and V?V,,,, are known numerically from the previous work. 
In the region inside S we have 


27/2 
[VV ,+-V°V,.—V3]= 5-H, . . . . (17) 


and we must solve these equations subject to the boundary conditions 
that 6 is finite at the impurity and the vacancy, tends to zero at infinity, 
and is continuous together with first derivatives on the surface S. 

As the potential has axial symmetry, the problem reduces to a two- 
dimensional relaxation solution. (For detailed discussion of two- 
dimensional relaxation methods, reference can be made to Southwell 
(1946) and Allen (1954).) Taking cylindrical polar coordinates with the 
impurity and vacancy on the z-axis and the origin halfway between 
(see fig. 3), the finite difference equations can be easily set up. The 
fact that the boundary is at infinity presents no real difficulty as 6 falls 
off quickly and may be put equal to zero to sufficient accuracy on a finite 
surface. Because of symmetry, we are only concerned in the relaxation 
solution with the half of the (p, z) plane on either side of the z-axis, the 
intersection of the finite surface with the (p, z) plane is then only part 
of the boundary ; the z-axis, which is the axis of symmetry, making up 
the other part. The complicating feature in the relaxation process is 
that we have a variable boundary S separating regions in which two 
different differential equations must be solved. The unknown boundary 
must again be found by trial and error, the boundary being adjusted at 
each stage of the calculation, so that eventually the function 6 is contin- 
uous with its first derivatives across S and the perturbing potential energy 
V, on the boundary is equal to the Fermi energy H,,. In practice, the 
condition (V,,),=#,, is used first to obtain a solution, then this solution 
is checked by the condition of continuity which is applied by considering 
derivatives along the mesh lines at the points where they cut the bound- 
ary S. Fortunately, however, the change in the boundary when the 
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impurity and vacancy are brought together is rather small in the present 
case. 
§ 6. Lyreraction Enercy Between IMPURITY AND VACANCY 


The interaction energy between the impurity and the vacancy in the 
metal is obtained by computing the difference between the energy of the 
impurity-vacancy pair separated by the assumed distance of 5 atomic 
units and the energy when they are an infinite distance apart. We shall 
consider separately the following contributions to the energy : (i) kinetic 
energy ; (ii) potential energy; (iii) exchange energy; (iv) correlation 
energy. 

It is certainly true that the method we have used to compute the 
potential neglects exchange and correlation effects. Thus the calculation 
neglecting exchange and correlation energies is more internally consistent. 
However, due to the variation principle of the Thomas—Fermi approxi- 
mation we should be able to calculate the energy including exchange 
to good accuracy by using the unmodified density to calculate the exchange 
energy. This is the procedure we shall adopt here, and we shall compute 
the correlation energy by a similar method. This latter term, in fact, 
turns. out to be relatively unimportant, but the exchange term is by no 
means so. 

The following expressions result for the changes in the various energy 
terms : 


(i) Kinetic energy change : 
5 2 ; 
AT = = (302)? | [n53+4+n 58,53 —n,5/3]dr  . . (18) 


where n, n,; and n, are the electron densities corresponding to the 
impurity—vacancy, isolated impurity and isolated vacancy cases res- 
pectively. 


(ii) Potential energy change : 


hoe eel 
AV=— 3 is [n—n,] dr+ 3 = [n—n,] dr 


1 : 1 
Ab Dy i [ V p(%— Ng) — V ps(M;—Mq) — V yy(My—Np) | dr— r. : =a (19) 
0 
(ui) Exchange energy change : 


oa we 
44=—3(=) | [n*3+-ng*/?§—nA8—n,¥3) dr. . . (20) 


eee a energy change 4W: (see especially Gombas 1949, 
ines 1955). 


This can be obtained in a quite straightforward way using the 
approximate formula 
0-0564n4/3 
niSLOOTR oe een 


for the correlation energy density as a function of the electron density n. 
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The change in energy in each case was evaluated by numerical integra- 
tion from our relaxation solutions for the densities and the values obtained 
were as follows: 

AT=—0-0042 
AV=+0-0011 | 
AA=+0-0018 
AW=-+0-0001 


(22) 


The total change is thus —0-0012, and so the energy is lowered by 
0-03(3) ev with respect to the separated imperfections. If we neglect 
exchange and correlation, however, the binding energy is increased very 
significantly, the value being 0-08(4)ev. It should, of course, be 
remembered that we have only included exchange in an approximate 
way even within the present framework ; it would be preferable to use 
the Thomas—Fermi—Dirac method from the beginning and calculate the 
total energy in this approximation. With the experience gained in the 
present investigation this should be possible, although it will undoubtedly 
be a very laborious task. We do not expect the interaction energy to 
change greatly from our result, although of course it should be remem- 
bered that the separate energy terms may well change appreciably. 
Correlation, as we have seen, is very much less important than exchange, 
and can probably be safely ignored in future work. 


§ 7. CONNECTION WITH FIRST-ORDER RESULTS 


As we have already discussed, Lazarus uses the simple form for the 
interaction energy given by eqn. (2) in his theory, with considerable 
success, and it is therefore of interest to compare the results of our 
computations with the first-order arguments of §3 which were shown 
to lead to eqn. (2) in the impurity—vacancy case. Putting Z=1, Z’=—1 
and 7)=5 in eqns. (7) and (8) and evaluating the kinetic energy change 


q° 
re ive V3 dr 
in confocal elliptic coordinates as 


Zexp (—a") me ey SE Gee el eo) 


the following results are obtained 


AT =—0-0040| (24) 
AV=+0-0023f ’ ae , 


to be compared with the results given in eqn. (22). 

It will be seen that both methods agree that there is a lowering of the 
kinetic energy which is sufficiently large to outweigh the increased 
potential energy, although the binding energy obtained from our compu- 
tations without exchange and correlation. (0-08(4) ev), is considerably 
larger than that given by the first-order method (0-045 ev). Further, 
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it is clear from our numerical results that the introduction of exchange 
markedly reduces the binding energy ; a conclusion which can again be 
qualitatively supported by a first-order argument. It is also of some 
interest to remark that the final numerical value we obtain for the 
binding energy (0-03(3) ev) is very close to that given by eqn. (3): 
(0-035 ev). 
§ 8. DIscuUSSION 
(i) Field round Impurity—-Vacancy System 

It seems of interest to note briefly here the main conclusions concerning 
the field and the distribution of electrons which can be drawn from our 
numerical calculations. We find that a simple superposition of the exact 
Thomas—Fermi potentials obtained for the isolated imperfections is a 
fairly good approximation to the potential for the two-centre problem. 
Also, as we have discussed earlier, using this model together with the 
Thomas—Fermi approximation there is a slightly distorted sphere round 
the vacancy (a little smaller than that for an isolated vacancy) in which 
the electron density is zero, and the negative charge is thus being shielded 
to the maximum degree by the unneutralized positive charge. 


(ii) Possible Errors Introduced by the Thomas—Fermi Approximation 

It should, of course, be borne in mind that the Thomas—Fermi 
approximation is not good in regions where the potential is varying 
rapidly, and corrections ought really to be applied. Attempts to do so 
have been made previously in simple cases (Scott 1952, March and 
Plaskett 1956) but so far no successful method having any degree of 
generality has been proposed. However, we feel that in the present 
problem, it might eventually be worth while to examine the Weizsacker 
correction to the Thomas—Fermi methody. 

It is, of course, difficult at the present time to assess the extent to 
which the errors in the procedure adopted here will affect the interaction 
energy. In this connection, however, we might perhaps remark that it 
is difficult to see how the present method applied in an unmodified form 
to two vacancies, could lead to results such as those of Seeger and Bross 
(1956) to which we referred earlier. It seems to us that the failure of 
the present method to allow electrons into regions of negative kinetic 
energy may be very serious in this latter case, and the possibility that 
this is a serious source of error in the present work cannot therefore be 
entirely excluded. 

§ 9. CONCLUSION 

We think, however, that in spite of the over-simplified model used 
and the limitations of the Thomas—Fermi approximation, the present 
result should be taken as giving some support to the approximate 
expression (2), perhaps with the modification of (3) introduced, as a 


+ Since this was written Dr. J. Friedel has informed us that he has carried 
out some computations of the Weizsicker correction (to be published). 
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useful representation of the interaction energy between an impurity and 
a vacancy in a monovalent metal. Thus, while further work is clearly 
necessary before a definite conclusion can be reached, we suggest that the 
present investigation tends to substantiate further the Lazarus theory 
of solute diffusion in metals. 
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ABSTRACT 


The graphical determination of half-lives of radioactive isotopes can be 
made more accurate if the remaining number of counts, instead of the 
counting rate, is plotted as a function of time. This method has been 
used to measure the half-life of 2’Si, for which the value 4-14-+-0-03 seconds 
has been obtained. 


§ 1. INTRODUCTION 

THE decay of a radioactive isotope is usually followed by counting the 
number of disintegrations (or a constant fraction thereof, according to 
the geometry and sensitivity of the counter) during successive time 
intervals. From the data so obtained the half-life is estimated either 
numerically or, less accurately, graphically from a semilogarithmic plot 
of the counting rate. It will be shown that the accuracy of the graphical 
method can be improved if a modified way of plotting is used. For 
comparison some numerical methods will be reviewed first. 


§ 2. NumericaL MetHops 
Peierls (1935) has given the following formula for the mean life time 
(7=t,/,/In 2) of a simple activity, if no background is present : 
7=t+t,[exp (t,/7)—1}-}, SS. Se eee 
where ¢, is the total counting time and # the time coordinate of the centroid 


of the observed counts. If 0=t)<t,<...<#, and if m, is the number of 
counts during the time interval 4t,=t¢,—t,_,, thent = >*n,t,/n, where 


n = >in, and#é, is the average time coordinate of the centroid of the 
counts in the interval J1,, i.e. 


t,=4(,-1+4,)—7[($.4t,/7) coth ($4¢,/7)—1] 
i 
=H(t,_+t,)— At, E (At, [r)— “5 Cine ee | 


+ Communicated by the Authors. 
{ Present address: Arthur D. Little, Inc., Cambridge, Mass., U.S.A. 
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Another formula, less simple than (1), is obtained if the method of 
least squares is applied to the semilogarithmic plot of the counting rate 
(fig. 1). The result is (Behrens 1951) 


DRL (t* —t, *)2 
ee a2 SUD aa ae (2) 
>kn,(t* —t,*) In (n,/At,) 
where ¢* = }‘n,t,*/n, and t,* is the time at which, in the average, the 
counting rate is equal to n,/At,, i.e 
t,*=4$(t,_1+4,)—7 In [sinh (4 4, /7)/( me 


1 


The calculation of the corrected mid-point ¢,* can be simplified if the 
time #, 4/2 which halves the number of counts, m,, in the interval Jt, is 
also recorded. The point ¢,* lies, on the average, very nearly two-thirds 
of the way from the mid-count point to the mid-time point so that 


t*=t, yet S341 ty) ty_s2ol=3(t,a1t+t_1/2+4,). cians (3) 
(The next term of this approximation is — 4t,/7207°.) This method has 
been found useful for plotting counting rates. 


Big. 1 


COUNTING RATE (CPS) 


10 10 20 
TIME (SECONDS) 


The mean square error c,?, due to statistical fluctuations, is the same 
for both expressions, (1) and (2), if 4t,<r. If moreover ¢, > 7, then 


(ae Ol ae ee ee (4) 


1000 I. Kuséer, M. V. Mihailovié and E. C. Park on a 


This is the minimum value of the mean square error which can be attained 
with the same total number of observed counts (7) (Peierls 1935). 

Instead of plotting the counting rate we may plot its integral, i.e. the 
number of remaining counts which at t=t, is N,= >*tin, (fig. 2). For 
this purpose the counting must be continued until the activity dies out or 
until the remainder (V,—7,,;=N )—n) can be estimated. The method 
of least squares now leads to 

ka 
dia, In (NQ/N,) - 

The weights a, must be so chosen that the mean square error becomes 
least. If At,<z, the proper choice is a,=4(n,+n,,,) for v<k, and 
A,=3N,+N441, Which minimizes the mean square error to o,2=77/n if 
t,>7. With such weights the formula (5) differs only insignificantly 
from (1). 


Fig. 2 
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§ 3. A GRAPHICAL METHOD 


As indicated above, the method of least squares, applied either to the 
numbers , or N,, leads to results of the same accuracy in both cases. 
The graphical method, however, distinctly favours the plot of the 
remaining number of counts (fig. 2), where it is easier to guess a correct 
straight line than in the counting-rate plot (fig. 1). Even a straight line 
through only two points, a few half-lives apart, will give good accuracy 
According to Peierls (1935) the points must be taken at t=0 and t= pen 
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in order to minimize the mean square error to o,2=1:5472/Ny. ‘This can 
be improved with three points (A, B, C) and a straight line which goes 
through the first point (A) and bisects the angle BAC. The minimum 
o,2—=1-2277/N, is reached if the points are taken at t=0, t=1-07 and 
t=2-57r. Hence any straight line, reasonably well fitting the points in 
the plot, will give the life time to an accuracy of the same order as the 
numerical methods, provided that the accuracy of the drawing itself is 
sufficient. 

Some details of this method will be included in the subsequent report 
about a measurement of the half-life of 27Si. 


§ 4. THE HALF-LIFE OF 2’Si 


This isotope was produced by the ?’Al (p,n) reaction in the external 
proton beam (10 Mey) of the Birmingham cyclotron, and later also by the 
8Si (y, n) reaction, using X-rays (maximum 26 Mev) from the Ljubljana 
betatron. In both cases a pneumatic tube system was used to transfer 
the target to a scintillation counter, in order to observe the decay. A 
chronograph, connected tothe output of a scaler, recorded every hundredth 
or thousandth pulse of the counter upon teledeltos paper tape (Kus¢cer and 
Mihailovié 1955). The scaling ratio was so chosen that more than a 
sufficient number of marks were obtained. When reading the marks a 
wider ratio (e.g. 1 : 200 or 1 : 2000) was chosen, so that in the central part 
of the diagram about two points per half-life were obtained. 


No =O 
(estimated) 


The further working scheme is shown by the accompanying table, which 
contains the data of a single experiment with a rather poor number of 
counts. The time coordinates ¢, of every second mark, i.e. of every 200th 


rd 
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pulse, are read from the paper tape (column 2). The differences At, 
(column 3) are taken, and the background and dead time corrections 
estimated (the sum of both being 6,, column 4). At the lower end of the 
table, where the background starts to exceed the true activity, the 
remaining number of ‘ true’ pulses (V,) has to be estimated. For this 
purpose an approximate value of the half-life can be used. Then the 
remaining counts N, (column 5) are built up, working backwards, by 
summing the n,=200—b,. The equation N)=*.200+N,— 4b, provides 
a check. Then the plot is made (fig.2) and the half-life read from it. 

Only a single decay of the 2’Al (p, n) activity, comprising over 90 000 
counts, has been used for the half-life determination (Park 1955). Some 
uncertainty was caused by an unknown dead time of the electronic 
equipment. The counting rate plot led to the value t,,,.=4:1-0:1 sec, 
and also permitted the dead time to be estimated. Using this estimate 
the data were replotted on a graph of the remaining number of counts 
(which is less sensitive to dead time variations), whereby the value 
t1.=4:14+0-03 sec resulted. 

Eight measurements with x-ray activated silicon were made, with the 
number of counts (Vy) ranging from 3000 to 7000, the total number of 
counts being 45000. The average happened to be identical with the 
previous result : 

t1/g=4 140-03 sec. 


The differences between the individual values of the half-life and the 
average were consistent with the expected statistical deviations. There- 
fore the standard error of the average was estimated according to eqn. (4), 
increased by the appropriate factor on account of the back-ground 
(Peierls 1935). The quoted result is in satisfactory agreement with the 
value 4-05-+0-10 of Hunt et al. (1954). 
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ABSTRACT 


It is established that the rectilinearity of the pits formed by etching 
octahedral faces of diamond depends in a sensitive manner on the tem- 
perature and rate of etch. Etching in potassium nitrate at 475°, 
although slow, produces rectilinear pits. At 525°c etching is much 
faster and the pits round off in the familiar manner already long known. 


§ 1. INTRODUCTION 


SEVERAL investigators have described the etching effects on diamond 
produced by strong oxidation. (For literature see Pandya and Tolansky, 
Proc. roy. Soc. A, 225, 33-48, 1954.) Thus etching can be produced by 
hot flames rich in oxygen, or by heating diamonds in a melt of potassium 
nitrate. In this laboratory extensive studies have been made with 
nitrate fluxes but at temperatures varying from 525°C to 750°c which 
are much lower than those used by earlier workers. 

On octahedral faces the etch pits are invariably equilateral, triangular 
shaped pits which are always oppositely oriented to the well known 
equilateral growth trigons. The pits so found are most frequently 
rounded off at the corners, whilst, on the contrary, the familiar growth 
trigons are almost invariably strictly rectilinear in outline. 

Now recently, in a personal communication, M. Omar of Cairo informed 
us that he and Kenawi have produced rectilinear etch pits by a novel 
method of etching{. Their diamond was supported on a tungsten 
filament at some 1000°c, in a vacuum chamber which contained a very 
small amount of air. After some hours, small, but rectilinear, etch pits 
formed, being, as always, oppositely oriented to the growth trigons. 

In this present paper we prove that the shape of diamond etch pits 
is determined by the rate of etching ; the pits can be arbitrarily produced 
either of rectilinear or curvilinear outline. 

ON OR ae DS eee 


+ Communicated by the Authors. 
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§ 2, EXPERIMENTAL 


We have carried out several experiments but select one set only. 
A diamond was cleaved on (111) in order to secure (111) surfaces free 
from the complication of natural trigons. We have already established 
(and are describing elsewhere) that when oppositely matched cleavage 
faces are etched, there is a remarkable correlation of the triangular 
etch patterns on the two separated faces. 

In etching with hot nitrate in this laboratory, we have usually formerly 
found that etching only really effectively sets in after 525°c, and have 
generally found that a one or two hour etch at this temperature produces 
a satisfactory early stage etch pattern, despite the fact that earlier 
authors advocated temperatures of at least 900°C and more. Etching 
below 525°c is normally so very slow that we have generally ignored 
lower temperatures in most of our etching experiments. However we 
began to etch here at 475°c. 

Etching was very slow indeed but we found that at this temperature 
the pits formed are in fact rectilinear (although small). Thus Plate 37 (a) 
(712) shows a region etched for no less than 24 hours at 475°c. Plate 
37 (b) (X712) shows the same region after being subject to attack for 
42 hours at 475°c. Plate 37 (c) is a selected region of this, at magnification 
1120. We see quite clearly that rectilinear pits have formed, increasing 
in size with time of etch. There is merely the barest suspicion of rounding 
at the corners in some of the pits, but not in all. 

Now a mere increase in etching temperature of only 50°c (to 525°) 
makes a striking change. We have traced this effect by two independent 
methods. In the first method we have etched the one face for a long 
time at 475°c. Simultaneously we have etched the other (untouched) 
matched cleavage face at 525°c for successive (much shorter) times. 
In the second method we have selected some of the rectilinear pits already 
formed at 475°c and then exposed them to the etchant at 525°. 

The slowly etched face of method (1) is that already shown in 
Plate 37 (a, b, c). The corresponding matched face is shown in Plate 37 
(d, e, f). Plate 37 (d) (x 712) shows that 33 hours at 525° leads to the 
long familiar pits which are distinctly rounding off into hexagon shapes. 
Plate 37 (e) (712) after 6 hours at 525°c shows the increase in size and 
in Plate 37 (f) at x 1120 a shape which may now well be compared with 
the rectilinear pattern of Plate 37 (c) which is shown by the matched 
cleavage etched at the lower temperature of 475°c. (It is of course 
quite clear that because we are dealing with a cleavage pair, the etch 
pits must of necessity be oppositely oriented on the two oppositely 
cleaved faces.) 

Equally striking are the effects shown by method (2). Thus Plate 38 (a) 
(x 712) shows a rectilinear pattern produced by a 42 hours etch at 
475°c. Plate 38 (b) (x712) then shows this identical region after a 
further exposure to etchant for 14 hours but now at 525°c. The marked 
set-in of rounding is clear. Further etching for another 13 hours at 
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525°0 leads to Plate 38 (c) (x 712) and here the long familiar rounded 
etch pits with flat bottoms have finally appeared. The growing pits are 
now beginning to eat into each other. This change in character has been 
produced by a mere 50°c increase in temperature. 

Another example of the successive conversion of a rectilinear system 
into a curvilinear system is shown in Plate 38 (d, e, f). Both examples 
in Plate 38 are typical of many cases we have studied. 

It should be recognized that even at 525°c the etching is still quite mild, 
for it still takes no less than 6 hours to produce the quite small pits 
of Plate 37 (e). Indeed it was this very slowness even at 525°c which led 
us initially to adopt this temperature as a reasonable starting place for 
our earlier extensive work on etching numerous diamond crystals. 


§ 3. CONCLUSIONS 


We have thus established two features (a) the rectilinearity or curvi- 
linearity of the etch pits depends only upon the temperature of etch, 
i.e. primarily on the rate of etch (b) rectilinearity is highly sensitive to 
temperature, a mere 50°c change altering the whole character. 

It is now clear that the appearance of rectilinear patterns merely 
implies a slow rate of etch. The rectilinearity of slowly formed etch pits 
has of course no connection at all with the familiar growth trigons which 
are oppositely oriented. 

In previous communications Tolansky only obtained the familiar 
somewhat rounded etch pits. We see now that this was because his 
etching, even at 525°c, was not sufficiently mild. It requires a really 
slow etching to retain rectilinearity. How slow an etching is needed is 
shown by the fact that the long 42 hour etch at 475°c produces pits whose 
sides are a mere 1/2000 cm in length. 
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§ 1. INTRODUCTION 


We have recently reported measurements on the y-rays emitted by the 
1-, T=1 state at 8-06 mev in !4N (Broude et al. 1957). This level has a 
rather large isotopic spin impurity (~7%) and a possible explanation 
of this has been given by Wilkinson and Clegg (1956) in terms of mutual 
contamination by the four states formed from a single 2s particle with the 
ground state of 130 as a unique parent. In a previous paper we suggested 
a possible identification of these four states on the basis of the weak 
decay modes of the 8-06 Mev state. The 8-70 Mev state in “N is, like 
the 8-06 Mev state, a very broad level formed by s-wave capture of 
protons by #8C. It is a 0-, 7'=1 state and possibly a member of this 
quartet of states. The present investigation had two purposes. Firstly 
to investigate the y-ray decay scheme of the 8-70 Mev level and secondly 
to investigate the isotropy of the ground state transition from the 8-06 Mev 
state to see if the ingoing particle is purely s-wave. If there is a substan- 
tial d-wave component present the number of similar states likely to 
mix in is increased and the whole situation is much more complex. 


§ 2, EXPERIMENTAL 


The 8-70 Mev state in N is formed by bombarding !8C with protons 
of 1:25 Mev. This energy is above the range of the 1 Mev H.T. set in 
this laboratory, but due to the large width of this resonance the yield 
at 1 Mev is one half the yield at the peak. Owing to this large width 
it is impossible to eliminate contamination y-rays by bombarding the 
target on and off resonance. Consequently we observed the y-rays 
at two different bombarding energies, 1 Mev and 900 kev, and compared 
the relative intensities. In addition the target backing was bombarded 
after each run as a further aid to the removal of background. A target 
of 70% C was used, of about 50 kev thickness at 550 kev and 30 kev 
thickness at 1 Mev. About 0-8% of the target area was found to be 
thick, giving rise to a small yield from the 550 kev resonance at 1 Mev 
bombarding energy. This contribution could be readily subtracted. 

The y-ray intensity measurements were carried out in the manner 
described in the previous paper. The relative intensities of the y-rays 
observed at the two bombarding energies after correction for the yield 
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from the thick target area are shown in table 1. The suggested decay 
scheme is shown in the figure. It will be seen that there is good agreement 
between the intensities at the two bombarding energies except for the 
y-ray leading to the 3-95 Mev level and its subsequent decay to the 
2-31 Mev level. As the 2-31 Mev level is fed in other ways the discrepancy 
does not show up so strongly in the intensity of this y-ray. 

Angular distribution measurements were made at three bombarding 
energies, 500 kev, 572 kev and 672 kev corresponding to mean energies 
in the target of 480 kev, 550 kev and 650 kev. The results are shown 
in table 2. 


Table 1 
Ep=1 Mev Ep=900 kev 
Relative Relative 
ayaa) intensity ste) intensity 

8-46 100 + 5 8-39 100 + 5 
5-69 3-2 + 0-4 5:69 2-9 + 0-4 
5-10 2-1 + 0-4 5-10 2-8 + 0-4 
4-92 3:2 + 0-4 4-9] 4-0 + 0-4 
4-5] 2-1 + 0-4 4-43 6-1 + 0-4 
3°55 4140-4 3°48 4-5 + 0-4 
3°38 6-5 + 1-0 3°39 8-9 + 1-0 
2-79 12-3 + 1-0 2-72 14-0 + 1-0 
2-31 79+1-0 2-31 10-5 + 1-0 
1-64 1-9 + 0-3 1-64 4-6 + 0:5 

Table 2 


Observed distribution 


Proton energy 


480 kev 1-+(0-00-£0-01)P, 
550 kev 1+ (0-008 --0-008)P,, 
650 kev 1+ (0-023 +.0-006)P,, 


§ 3. Discussion 


In addition to the main ground state transition, direct transitions 
from the resonance level occur to the levels at 3-95 Mev, 4:91 Mev, 5:10 Mev 
and 5:69 mev. These are the same levels as are fed from the 8-06 Mev 
level, with the exception of the 5-10 mev level which was not fed from 
the lower resonance. 

The strongest of these transitions is to the 5-69 Mev state. In the 
previous paper we suggested that this state might be the 1”, 7—0 
member of the quartet of s-particle states, as its decay to the ground and . 
first excited states is in the expected ratio if it is the sole contaminator 
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of the 8-06 Mev state. Since that work was completed Wilkinson and 
Bloom (1957) have shown that the 6-23 mev level in N also decays to 
the ground and first excited state in the correct ratio for it to be the sole 
contaminator for the 8-06 Mev level. This contaminating state should 
be 1- and the observed strength of the transition to the 5-69 Mev state 
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from the 1-25 Mev resonance level helps to decide whether this is a possible 
spin and parity for the 5-69 Mev state. Taking the radiation width of 
the 8-70 Mev level for decay to the ground state as 46 ev (Woodbury et 
al. 1953) the radiation width for decay to the 5-69 Mev level is 6-0 oa 
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This corresponds to 13 Weisskopf M1 units, which is four times higher 
than the strongest authenticated M1 transition (Wilkinson 1956). If 
this transition is an E1 transition then the radiation width is 0-45 Weiss- 
kopf units, which is above average El strength. Thus it seems fairly 
certain that the transition to the 5-69 Mev level is an E1 transition and 
this state cannot be 1-, but is 1+. 


3.1. The 4:91 MeV Level 


The existence of the transition to the 4-91 Mev state from the 0- 
8:70 Mev state shows that the former level cannot have zero spin, as is 
generally assumed. The evidence for this assignment comes from the 
8C (d, n)'4N stripping results of Benenson (1953) which appear to show 
that the 4-91 mev level is 0 or 1~. The absence of a transition to the 0+ 
state at 2-31 Mev and the fact that the 4-91 Mev level is not fed from the 0+ 
resonance at 1-16 Mev have previously favoured the 0~ assignment. 
The existence of the transition to the 4:91 Mev state at this resonance 
means that the level is presumably I~. However, the decay of this 
state is difficult to reconcile with the 1~ assignment, for the isotopically 
allowed E1 transition to the 0+, 2-31 Mev state is not observed whereas 
the isotopically forbidden transition to the ground stateis. The conditions 
for observing these transitions are better at the 550 kev resonance but 
here also the isotopic spin favoured transition was not observed. Further 
the radiation width for the transition to the 4-91 Mev state is 1-8 ev 
which corresponds to the rather high value 1-9 Weisskopf M1 units. 
It thus seems that the 4:91 Mev state requires further investigation 
before its spin and parity can be stated with confidence. 


3.2. The 5:1 MeV Level 


The radiation width for the transition to the 5-10 Mev state corresponds 
to 1-4 M1 units or 0-047 El units. These transition strengths are consis- 
tent with an assignment of 1~ or 1+ but the absence of a strong transition 
to the 2-31 Mev level from the 5-10 Mev level rather favours I+. 


§ 4. ANGULAR DISTRIBUTION 


The angular distribution measurements show some evidence for 
anisotropy at 550 kev bombarding energy and a quite definite anisotropy 
at 650 kev. At 650 kev the ground state yield is partially due to a 
contribution from the very wide 1-25 Mev resonance which is essentially 
isotropic and adds incoherently as the 1-25 Mev and 550 kev resonances 
are formed through different spin channels. In addition there is also 
the small thick target yield. After correcting for these effects we are 
left with a distribution of the form 1+-(0-04-+-0-01)P,. Using the radius 
which Milne (1954) has used for the analysis of the elastic scattering 
results, namely r—(1-+-1-41 13/3) x 10-8 em we have computed the 
amount of d-wave contribution required to produce such a distribution. 
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The relative phase of the ingoing s- and d-waves at these bombarding 
energies is determined almost exclusively by the coulomb field and is 
very nearly independent of the assumed value of the nuclear radius. 
At 650 kev this angle is 81° and is increasing slowly as the bombarding 
energy is lowered. As a result, the interference term, though still the 
dominant term, is small and the d-wave intensity required is 6-05% 
of the s-wave intensity. The relative intensities of the s- and d-waves 
calculated from barrier penetrability alone is 1-69%. 

On the assumption that the ratio of the s- and d-wave intensities 
varies across the resonance in the ratio of their penetrabilities, the expected 
distribution at 550 kev is P,+0-014P, in reasonable agreement with the 
observed 1+(0:008-+0-:008)P,. The reduced width of this resonance 
is 13°, of the Wigner s-wave limit. The d-wave contribution is thus 
approximately 45°/, of the Wigner d-wave limit. Any reduction in 
the radius used makes the d-wave contribution an even larger proportion 
of the d-wave sum. It thus seems that the description of this level in 
terms of the °C core and a single 2s particle is not well founded and that 
mixed configurations must be employed involving a large amount of 
d-wave. 
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ABSTRACT 
It is shown that there is a strong correlation between flow stress and 
subgrain size in polycrystalline aluminium deformed in tension; the 
strength of the boundaries does not appear to depend on boundary 
angle. The variation of flow stress with temperature suggests that the 
strength controlling factor is elastic interaction. 


§ 1. INTRODUCTION 


KELLY (1954) extending a theory of Gay et al. (1954), has suggested that 
the flow stress, o, at any point on a stress/strain curve of a polycrystalline 
metal specimen is determined by the instantaneous value of the subgrain 
size, t, according to the relation ooct-'/?. Using the experimentally 
observed variation of subgrain size with strain he finds good agreement 
with the stress/strain curves of aluminium and iron. 

However, determinations of subgrain size at low strains, where most 
change occurs, show considerable scatter, and the effect of other variables 
which might be expected to affect the flow stress, such as boundary angle, 
was not investigated. For these reasons the subject has been re-investi- 
gated. Measurements have been made, using the x-ray microbeam 
technique, of the subgrain size and range of misorientations in one original 
grain in a series of specimens, and these have been correlated with the 
flow stress. For details of the x-ray measurements see Hirsch and 
Kellar (1952). In addition, the variation of flow stress with temperature 
has been determined, in the hope that it would throw some light on the 
processes occurring at yield. 


§ 2. EXPERIMENTAL METHOD 


Tensile test specimens were made from superpure aluminium (99-994'/) 
and recrystallized at 550°c to give strain free grains of size ~200 x. 
These specimens were then strained by amounts ranging up to 35% 
extension at temperatures in the range —183°c to 375°c. By this 
means a set of specimens was prepared in which subgrain size and 
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misorientations were not closely related, as would have been the case for 
deformation at a single temperature. The influence of temperature of 
deformation on substructure will be described in a later paper. 


Fig. 1 


STRESS f 


Kgm/mm* 


STRAIN Timm 
Typical restrain curve. 
Before x-ray examination all specimens were etched to a depth of 


fe ft, as the surface layers are not characteristic of the interior of 
1e metal. Microbeam photographs showed that internal strains in 


Flow Stress of Polycrystalline Aluminium 1013 


specimens deformed at low temperatures were large. Accordingly, these 
specimens were annealed for 12 hours at 140°c. As Kelly (1953) has 
observed, the effect of such a treatment is to sharpen the spots on a 
microbeam photograph and to eliminate the background between them, 
without greatly affecting either the subgrain size or the angular extent 
of an arc. After this treatment internal strains, as estimated from 
microbeam photographs, were comparable in all specimens. 

The specimens were finally restrained at 0°c in a hard beam machine 
of Polyani type at a rate 2-5x10-> sec. Typical restrain curves are 
shown in fig. 1. The intersection of the two straight portions of each 
curve is taken as a measure of the flow stess, o. Some plastic flow 
does occur at lower stresses, particularly in specimens that have been 
annealed. In table 1 are listed values of subgrain size, ¢, range of 
misorientations, 8, and flow stress, o, for all specimens. 


Table 1 

tu B° o Kg/mm? tp B° o Kg/mm? 
40 1-0 1-04 3-6 (Fy 3:19 
25 1-0 1-04 4:3 4-5 3-53 
18 2:3 1-44 3-3 7-0 4-2] 
14 4:5 1-40 9-0 1-0 2:17 
20-5 0-8 1-41 6-6 2-0 3-40 
16-5 2-6 1-7] 4-45 3-6 4-02 
12 4-0 2-56 2-5 7:8 4-34 
16 1-0 E91 3-2 4-5 3°61 
11-5 3-0 2-25 1-8 9-0 5-20 

5-2 5-0 2-74 1:0 7:3 7:36 


The variation of flow stress with temperature cannot be deduced from 
these results as many of the specimens were annealed before restraining. 
Instead, following Cottrell and Stokes (1955), specimens were deformed 
5% at a high temperature, cooled to the reference temperature and 
restrained to just beyond the flow point, strained a further 5% at the 
high temperature and so on, giving several values for the ratio of the 
flow stresses at each pair of temperatures. The temperatures were 
obtained with the following baths: 


—177°c solid CO, in acetone 
0°c ice in water 
100°c oil bath controlled to --1°c. 
200°c to 500°C salt bath controlled to +-2°c. 


The time taken to change a bath was standardized at 15 minutes. 
Results are summarized in table 2, together with the ratio of the flow 


« 
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stresses after correction for the change of elastic constants with temper- 
ature (taken from Koster 1948). A high reference temperature was 
chosen for specimens strained at the highest temperatures as the initial 
rate of work hardening on restraining at 0°c was too great to allow an 
accurate determination of the flow stress to be made. These readings 
are scaled to a reference temperature of 0°c using the relation 


(o9/3)(03/04)=(29/01)- 


This relation has been tested by Cottrell and Stokes (1955). In fig. 2 
(c/o), after correction for the change in shear modulus, is plotted 
against temperature. 


Table 2 
L056) T Rep © o7/T Ret o7/o9 
®) IT 0-94 0:95 0-95 1-02 
100 0) 0-96 0-89 0-90 0:97 
200 O 0-89 0-90 0:85 0-97 
250 (0) 0-79 0-80 0:77 0-89 
300 0 0-68 0-65 0-65 0-78 
350 0) 0-54 0-56 0-55 0-68 
400 300 0-73 — == 0-62 
450 300 0-66 = — 0-60 
500 400 0-92 — = 0-64 
550 450 0:90 — — 0-61 
Fig. 2 
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1 600 


Ratio of flow stress op at temperature 7’ to flow stress Oy at 0°c. 
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§ 3. Discussion 


Consider dislocations piled up against a subgrain boundary. The 
shear stress at the head of a pile up of » dislocations under an applied 
shear stress 7 is nr. If the criterion for flow is that this stress shall 
exceed some boundary strength 7+,, we have 


CTs . 


Now 


where / is the length of the pile up, k is a factor of order unity, G the 
shear modulus, and 6 the Burgers vector of the dislocations (Cottrell 


1953). Hence 
Gb\ 1/2 
— (=) [2p 1/2, 
The tensile flow stress o will be proportional to 7. If the polycrystalline 
specimen is homogeneously stressed, for the most favourably oriented 
slip system we have tc. The boundary angle « may be expected 
to affect t,. The relation between « and the experimentally observed 
quantity 6 is not known with certainty. For a random distribution of 
orientations (Gay et al. 1953), 
a~B/3, 
and this will be assumed. If we equate the length of a pile up to the 
subgrain size t we get a relation between o, ¢t and f of the form 


o=At-¥f(B) 
where A isa constant. ‘To test for this a regression equation of the form 
log c=a,+b, log B-+b, log t 
has been fitted to the data, by standard methods, with the coefficients 
a,=0-898, b,=—0-024, b,=—0-558. The coefficients of correlation are 
119-3 —0-08, 715-2=—0-92. It appears, then, that there is a strong 
correlation between o and ¢t which is not significantly different from that 
proposed by Gay et al. namely o oct-!” but that o and f are virtually 
uncorrelated. In fig. 3 o is plotted against ¢1/. 
The slope of the line, with the assumptions mentioned, yields 
Tp~G/50. 
It is not at present possible to explain the strength of the boundaries. 
Gay et al. considered the stress a distance « ahead of a pile up, and 
took as their criterion of flow the formation of a dislocation loop by 
thermal fluctuations. According to Frank (1950) this can occur at 
room temperature if the stress is of the order of G/25. In view of the 
value for the boundary strength deduced here, it may seem surprising that 
they found such large stresses to occur. Part of the explanation is that 
they took t=c. Also, as their specimens were not annealed, their flow 
stresses were higher than for comparable specimens in the present work, 
The anneal also reduces the subgrain size slightly. 
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The stress distribution inside the polycrystalline specimen is not 
known, and it may be that stresses at the head of a pile up are as high 
as G/25. Most of the energy to form the new dislocation loop is provided 
by the applied stress, so dependence of the critical strength on temperature 
is slight, and would probably not be detected. 

In the neighbourhood of 0°c the change in flow stress with temperature 
is only slightly greater than the change in the elastic constants (fig. 2). 
This suggests that if the criterion for flow is breakdown of a boundary, 
the strength controlling factor is elastic interaction. In this case, a 
strong dependence of boundary strength on angle would be expected, and 
no such dependence has been observed. This, however, may be because 
the relation assumed between « and f is incorrect. The range of f is 
not large, and the range of » may very well be much less. 


Fig. 3 


-1/2) 


(u 
Tensile flow stress o as a function of (subgrain size t)-1/2, 


An alternative explanation is that the leading dislocation of a pile 
up combines with the dislocations of the boundary. The greater the 
density of dislocations in the boundary, the smaller the length of combina- 
tion of each with the leading dislocation of the pile up, giving a very much 
smaller dependence of strength on boundary angle than would be expected 
in the absence of combination. If the boundary contains dislocations 
of the correct slip system, the combinations may take the form of 
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Cottrell-Lomer locks. The mechanism of failure, however, would not be 
that considered by Stroh (1956), as the strength is not temperature 
dependent. 

In any case, a dependence of strength on boundary angle would be 
expected for very low angle boundaries. This is being investigated. 

The drop in the curve of fig. 2 for temperatures above 200°C requires 
some comment. 200°c is the temperature at which grain boundary slip 
begins to operate in aluminium (Ke 1947), and it seems probable that 
this is the explanation. Relaxation of shearing stresses across grain 
boundaries would lead to a stress redistribution, and it may be that the 
effective internal shearing stresses at the places where slip is occurring 
are just as high as at low temperatures, when allowance has been made 
for the change in shear modulus. 
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ABSTRACT 


The nature of the damage in neutron irradiated lithium fluoride has 
been examined by x-rays up to neutron doses at which the majority 
of the ions have been displaced from their lattice sites. The observations 
have been made by Debye-Scherrer powder photography, Laue photo- 
graphy and small-angle scattering. The results indicate that isolated 
defects predominate at low doses, but that these cluster to form large 
defects at doses higher than about 10!” neutrons/cm?. 

With the exception of some small-angle scattering all the x-ray effects 
anneal out in the range 400°-500° c. This residual scattering is inter- 
preted as that associated with gas bubbles. 


§$ 1. INTRODUCTION 


THE nature of the damage induced in lithium fluoride by reactor irradia- 
tion has been the subject of several recent x-ray investigations. The 
problem has doubtless owed its popularity to the simplicity of the un- 
irradiated structure and, more important, to the ease with which the 
structure can be damaged. 

Binder and Sturm (1954) showed that for a dose of 6 x 1016 neutrons/em? 
the lattice expansion as measured from the change of lattice parameter 
was equivalent to that determined from density measurements. This 
indicated that the expansion was due to the production of isolated 
Frenkel defects, an interpretation which was later supported by annealing 
experiments (Binder and Sturm 1955). Keating (1955) measured the 
integral breadths of the 001 reflections of crystals irradiated at rather 
higher doses, 2-1 and 7-510! n/em?. He observed considerable line 
broadening, which he interpreted in terms of a model of flat clusters 
of defects in the (001) plane separated by uniformly stressed material. 
Finally Mayer (1955) has reported measurements by Perio of the variation 
of lattice parameter with dose, and Brooks (1956) has reported some 
x-ray results of Guinier. Both these sets of results will be referred to 
below. 


In the present work new x-ray results are reported on lithium fluoride 
irradiated in the range 3x 1016 to 7x 1018 n/em?. 
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§ 2. IRRADIATION PROCEDURE 


The standard samples were plates, } mm thick, cleaved from larger 
single crystals obtained from Harshaw Chemical Co., U.S.A. These 
crystals were irradiated in the centre of Brpo, the graphite moderated 
reactor at Harwell, in a slow neutron flux of about 1012 n/em2/sec 
and at a temperature of 55-L15° o. 

The doses were measured by determining the y-activity of cobalt wire 
monitors placed with the samples. Owing to self-screening the mean 
dose received by a sample was less than the unscreened dose recorded 
by the monitor. Calculations based on the assumption of an infinite 
slab of }mm thickness showed that the standard sample received a 
mean dose 80% of that recorded by the monitor. This correction has 
been applied to the dose figures quoted below. Self-screening also 
produced a variation of +3°% in the dose throughout a given sample ; 
this variation is about half the error involved in estimating the relative 
doses. 

§ 3. X-Ray RESULTS 


3.1. Measurements of Lattice Parameter and Line Breadth 


Debye-Scherrer photographs were taken of powders, which were 
prepared from fragments cleaved from irradiated single crystals. The 
lattice parameter a was obtained from the line positions on the photo- 
graphs. Figure 1 shows the variation of a with dose. For small doses 
the change of a, 4a/a, is proportional to dose ; at higher doses @ increases 
to a maximum at 41017 n/em? and then falls again. These results 
are in good agreement with those obtained by Perio (see Mayer 1955). 


Fig. 1 
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_ Lattice parameter of LiF as a function of slow neutron dose. 


The integral breadths f,,,, of the powder lines were measured from 
line profiles obtained by microdensitometry. The instrumental broaden- 
ing was corrected by the Jones (1938) method, and it was assumed that 
the breadth of the unirradiated crystal was entirely instrumental. A 
plot of £ cos 6/A versus 2 sin 6/A for the most heavily irradiated sample 
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(7-2 x 1018 n/em2) was a straight line of slope 0:0042 through the origin. 
This showed that the line broadening was predominantly due to strain 
rather than crystallite size, and that the mean strain was 0-0042. 
The integral breadth of the 422c, line as a function of dose is shown in 
fig. 2. B4o9 is the same as the unirradiated value up to Uso n/em? ; 
above this Bag increases rapidly with dose and shows no sign of saturation. 


Fig. 2 
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Integral breadth of 422«, line ©, and intensity of x-ray scattering at 1° x, 
as a function of slow neutron dose. 


3.2. Laue Photographs 


Transmission ‘Laue*® photographs, using both white and mono- 
chromatic MoK« radiation, have been taken of a number of irradiated 
single crystals. In the heavily irradiated (i.e. greater than 1018 n/em?) 
samples narrow, well-defined streaks appear along the <100) directions. 
The streaks extend from the origin to the 200 reciprocal lattice points, 
and are produced by reflection of the characteristic MoK« radiation. 

‘Defect diffuse scattering ’ (Huang 1947), which arises from the dis- 
tortion of the lattice around the individual point defects, has been carefully 
looked for. This scattering is difficult to observe experimentally, as it is 
concentrated near the strong Bragg reflections and is in a region where 
the thermal diffuse scattering is a maximum. However, extra diffuse 
scattering has been observed along the radial direction joining the 
reciprocal lattice point to the origin. Such a radial distribution of 


intensity is predicted by the theory of defect diffuse scattering (Huang 
1947, Cochran 1956), 
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3.3. Small-angle X-Ray scattering 


The apparatus used for these tests was that described by Smallman 
and Westmacott (1957) for the study of cold work and irradiation effects 
in metals. It consisted essentially of an asymmetric quartz mono- 
chromator limited with slits to give }° beam divergence. The sample to 
focus distance was 80 mm and the intensity of scattering recorded either 
on a film or directly using a Geiger counter. 

A film taken with a crystal irradiated to 3x 1018 n/em? and oriented 
with the (100) and (001) planes respectively perpendicular and parallel 
to the x-ray beam, showed that the small angle scattering was approxi- 
mately isotropic. It was therefore possible to use a Geiger counter 
(Mullard type MX118 with slit 0-3mm x5 mm) to measure accurately 
the intensity distribution along the (100)-(001) zone. Figure 3 shows 
a plot of log (intensity) versus 6? (radians?) for six crystals irradiated 
with successively increasing doses of 5x 101% 1:01017, 1:0 1018, 
3°3X 1018 and 7x101%n/cem?. The results (corrected for parasitic 
scattering but not for beam height) have been plotted in this way since 
it is known (Guinier and Fournet 1955) that for a collection of identical 
scattering groups the intensity scattered over an angle @ is given by 

I=Mn?I, exp (—47?K767/3A2). Rarer nur ners eh 
Here /, is the intensity scattered by an electron at small angles, M is the 
number of groups, ” represents the difference between the number of 
electrons in one group and in an equal volume of surrounding matrix, 
and K is the radius of gyration of the groups. 

The amount of small angle scattering for irradiations less than 
31017 n/cm? is extremely small but above this dose the scattering 
becomes appreciable. The intensity of scattering at 1° as a function 
of dose is shown in fig. 2 together with the line broadening results. 

Accurate analysis of such a complex scattering system is difficult but 
an estimate of the maximum size of scattering region is of interest. 
The table shows the radius of gyration K for the various irradiated 
erystals calculated from the extreme low angle linear portion of the 
log J versus 0? curve. 


Trradiation dose Size K(A) Irradiation dose Size K(A) 
(njem?) (n/em?) 
‘spo Mi be 10 Loc io. 2 
1x10" 15-5 a3 LOR 22 
3:7 10" 18-5 (ie WU 22 


It can be seen that the size of the scattering group increases with 
neutron dose up to about 10!8n/em? and then remains approximately 
constant at 224. The number of such regions also increases extremely 
rapidly with dose, as the intercept on the log J axis indicates, 
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Fig. 3 
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Log (intensity) versus 6? for six crystals irradiated at successively increasing 
doses. 


3.4. Annealing Experiments 


Irradiated single crystals have been vacuum annealed for 1 hr periods 
at 50° and 100° c intervals in the range from room temperature to 700° c. 
After each anneal the crystals were examined by the three x-ray techniques 
mentioned above. The bulk of the damage annealed out between 400° 
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and 500° c., ie. the lattice parameter reverted to its original value, the 
powder lines became sharp and the Laue streaks disappeared. The 
small angle scattering results showed that the annealing behaviour of 
the scattering group was extremely complex. Below 400°c there is a 
small increase in small angle scattering but a large decrease does occur 


Fig. 4 
10° 
| 400°C 
Qs irradiated 
i \\4 
ie | 
c lo rf 
= 
— 
2 
g \ 
Eile 
3 A \ 
~ \ 4 
£ ~B 
350°C 
10° See 
r? x E a 
r 500°G 
600°C 
10! i ! 


0 0.2 0.4 0.6 0.8 ne 1.2 
Q* x 10° —> 
Log (intensity) versus 6? for crystal irradiated to 1018 n/em? and annealed for 


one hour at successively increasing temperatures. @ As irradiated, 
[] annealed at 350° c, x 400° 0, A 500° c, gy 600° c. 


in the range 400-500°c. Further annealing at temperatures above 
500° c causes the log J versus 6? plot to become approximately linear but 
does not remove the scattering completely. The scattering curves 
for a crystal irradiated to 1x 10'8 n/em?* after various anneals are given 
in fig. 4. Substantial small angle scattering remains after high tempera- 


ture annealing. 
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Irradiated lithium fluoride is dark red or black depending on the severity 
of the dose. On annealing the colour of all the samples changed. Crystals 
irradiated less than 1018 n/em? become yellow-brown at low temperatures 
but finally colourless in the range 400-500° c. The crystals more heavily 
irradiated remained dark up to 400-500° c, then became white. 


§ 4. DIscUSSION 


In the irradiation range below 101? n/cm? the amount of small angle 
scattering and line broadening is extremely small whereas the lattice 
parameter increases linearly with dose. This initial linear portion of the 
lattice parameter versus dose curve (fig. 1) can be interpreted as being 
mainly due to isolated Frenkel defects distributed at random throughout 
the lattice. From the slope of this portion we can obtain, by assuming 
that the lattice expands by the molecular volume for each pair of Frenkel 
defects produced, an estimate of the number of defects per fission. This 
number is fairly close to the theoretical estimate of 1900 obtained by 
assuming that all the defects remain in the lattice during and after the 
irradiation (Seitz and Koehler 1956). We conclude therefore that the 
major part of the damage does not anneal out at room temperature. 

We can interpret the saturation of the lattice parameter versus dose 
curve by using the rough criterion for saturation suggested by Pease 
(see Varley 1955). An upper limit to the concentration of Frenkel 
defects is that for which any additional interstitial introduced in the 
lattice is adjacent to a vacant lattice site. If the defects are immobile 
this will occur when the vacancy concentration is of the order of 1/N, 
where VV is the number of lattice sites surrounding an interstitial site. 
Now the cross section for thermal neutron absorption in LiF is 70 barns 
so that for the saturation dose of 3x101!7n/em? a mole fraction 
(70 x 10-*4) (3x 1017)=2-1x 10° of fissions occurs. The corresponding 
concentration of vacancies is (2-1 x 10~5) (4) (1900)=1/50. This is 
rather less than 1/N, but the difference can be accounted for by assuming 
interstitial-vacancy interactions over a few ionic distances. 

We attribute the fall of the lattice parameter versus neutron dose curve 
to the coagulation of point defects. Evidence for coagulation is provided 
by the appearance of (100) streaks in Laue photographs and by small 
angle scattering. The streaks indicate a break-up of the periodicity 
in the (100) directions, as produced, for example, by the condensation 
of interstitial atoms in (100) planes. The small angle scatter curves 
show that the scattering system is more complex, being more sensitive to 
vacancy groups than interstitial aggregates. However it is expected 
that aggregates of lithium and fluorine atoms, together with groups of 
vacancies are present. Any variation in size of each of these scattering 
groups will give rise to a concave type of log (intensity) versus ‘62 curve, 
as observed. Coagulation of defects also produces long range strain, 
as revealed by the line broadening measurements. The interpretation 


X-Ray Study of Neutron Irradiated Lithium Fluoride 1025 


of line broadening in terms of crystallite size is invalid (see § 3), although 
Varley (1956) reports that Pease has deduced thereby an estimate of the 
size of the region disordered by a single fission process. 

On annealing, the Laue streaks, line broadening and changes in lattice 
parameter all completely disappear in the temperature range 400-500° c. 
However, some small angle scattering still remains even at 700° c. This 
residual scattering is believed to be due to the new atoms, 3H and ‘He, 
formed by transmutation. The number of such atoms per group can be 
calculated from the x-ray measurements, since from eqn. (1) the number 
of scattering particles M is given by 


M=NI(0)/n?, ays eee een a} 


where /(Q) is the intercept intensity in electrons per atom and N is the 
number of atoms/cm*. The conversion of the intensity from counts/sec. 
to electrons/atom can be made, knowing the geometry of the apparatus 
and the main beam intensity (Smallman and Westmacott 1957). For 
our arrangement one count/sec. is equivalent to 20 electrons/atom. 
Thus in the 1x10! n/cm? and high temperature annealed sample 
the number of groups is about 4 1015/em’. However, the number of 
transmuted atoms formed during irradiation is 8 x 101%/cm? so that the 
number of atoms per group is about 2000. If these small atoms fit 
interstitially in the lattice sites ($,4, 3) or ($,0,0) then the volume 
occupied is about 500 unit cells. This is not in agreement with the 
observed size of the scattering group of 304A radius. The difference 
can be accounted for if it is assumed that gas atoms by capturing vacan- 
cies in their neighbourhood expand into gas bubbles forming a stable 
cavity. Such behaviour has been observed metallographically when 
irradiated uranium (Cottrell 1956) and beryllium (Barnes and Redding, 
private communication) is heated. An estimate of the size of cavity 
which this group of atoms would acquire at temperature 7°K may be 
obtained since the pressure of a gas bubble of radius # is given by 


Pe meee ee ee (3) 


where y is the surface energy of the lithium fluoride. Since the volume 
of gas atoms at n.t.p. is 
eee ae en oe ees, (A) 


where p, is the density of the gas at n.t.p. and A, the atomic weight of the 
gas and N Avrogadro’s number, the number of gas atoms per bubble 
at pressure P and temperature 7K is given by 


Dees AGT 5 5° =~. «i(B) 


Taking the number of atoms per bubble as that given above (2x 10°) and 
with y= 200 ergs/cm? and pg/Aq=4:5 x 10-* g/cm? the size of the bubble 
is about 30 radius, in agreement with that observed at 600° c. 
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§ 5. CONCLUSIONS 


Three principal x-ray techniques, Debye—-Scherrer powder photography. 
Laue photography, small angle scattering, have been used to determine 
the effect of neutron irradiation on crystals of lithium fluoride. The 
results indicate that for low doses (<101!’ n/cm?) the defects are pre- 
dominantly isolated. However at doses higher than about 10!” n/cm? 
clusters of defects such as interstial platelets in (100) planes and groups 
of vacancies are formed. These defects anneal out in the temperature 
range 400°-500°c but above this temperature it is believed that the 
transmuted atoms 3H and *He which remain form small gas bubbles. 
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Note added in proof.—The presence of gas pockets has been confirmed 
by microscopic examination of samples annealed for long periods just 
below the melting point. For example, bubbles six microns in size and 
of square cross section were observed after heating for 24 hours at 700° c. 
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SUMMARY 

Mason’s theory of Rochelle salt is critically examined in relation to 
Bancroft’s high-pressure measurements of the dielectric constant. It is 
shown that, though the theory as it stands cannot explain the results, 
it can be made to do so fairly well by a slight modification. The theory, 
however, requires the effective dipole moment to vary with volume in 
an unexpected way, and no satisfactory explanation has been found for 
this. 


§ 1. [yTRODUCTION 

ROCHELLE SALT was the substance that was first shown to be ferroelectric. 
In recent years many other ferroelectrics have been discovered, but 
Rochelle salt has properties not possessed by any of the others. The 
most striking of these properties is the existence of two Curie temperatures. 
It is true that many other ferroelectrics have two or more transition 
temperatures, but each transition is always to a new crystal symmetry ; 
it is only in Rochelle salt that the symmetry? below the lower transition 
temperature is the same as above the upper transition temperature. 
This property alone would make it of special interest to develop a theory 
for Rochelle salt, but in addition we are able to check a theory more 
closely for Rochelle salt than for any other ferroelectric. This is for the 
following reason. Measurements of the dielectric constant have been 
made by Bancroft (1938) for various pressures up to 10 000 atmospheres. 
These measurements not only give us the dielectric constant as a function 
of pressure and temperature, and hence of specific volume and tempera- 
ture, but also give the transition temperatures as functions of pressure 
or specific volume. They show, incidentally, that there are two transition 
temperatures under conditions of constant volume as well as constant 
pressure, and also that if the substance is kept at constant temperature 
and the pressure varied there may be two transition pressures. 

The dielectric properties of Rochelle salt have been very fully discussed 
in a phenomenological way by Mueller (1940), but he did not attempt 
to relate them to the atomic structure except in so far as it determined 
the crystal symmetry. The most successful atomic model is that suggested 


+ Communicated by the Author. 
+ This does not necessarily imply that the structure of the unit cell is exactly 


the same. 
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by Mason (1947, 1950). He was able to account for the existence of two 
Curie temperatures at normal pressure, but did not investigate whether 
the model gave two Curie temperatures at constant volume as Bancroft’s 
results show that it should. We have therefore examined the model 
in some detail; we find that, in the exact form that Mason gave it, it 
does not give two Curie temperatures at constant volume, but that by 
a slight modification can be made to do so. With this modification we 
find that we can explain Bancroft’s results fairly well, but that some 
difficulties still remain. 


§ 2. EXPERIMENTAL OBSERVATIONS 


If the dielectric constant of Rochelle salt is measured at normal pressure 
as a function of temperature it shows two sharp maxima at about —16°c 
and 24°c. In between these temperatures the substance shows hysteresis 
and other properties characteristic of a spontaneously polarized body. 
The spontaneous polarization has been measured ; its maximum value 
is about 2-5 x 10-7 coul/em?, and it tends to zero at the two transition 
temperatures. The transitions are accompanied by changes in crystal 


Fig. 1 
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The reciprocal susceptibility of Rochelle salt as a function of pressure at 40°C 
(Bancroft 1938). 


form ; when spontaneously polarized the crystal has monoclinic symmetry 
(class Cy), and otherwise is orthorhombic (class D,). The transitions 
appear to be second order ones. Anomalies also occur in some of the 
piezoelectric constants at the transitions, but it has been shown by 
Mueller that these are the natural consequences of the anomalies in the 
dielectric constants. 

Bancroft (1938) has made measurements of the dielectric constant at 
various pressures up to 10 000 atmospheres and for temperatures from 
—20°c to 60°C. In fig. 1 we show his values of the reciprocal susceptibility 
plotted as a function of pressure at 40°c. The reciprocal susceptibility 
at constant stress along the a-axis (the axis of spontaneous polarization 
in Rochelle salt) is defined by the equation 


x_ (0h, 
NS SANGY EI). wh eine ee ae ae) 
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where X denotes the stress. At a second-order transition y,,* will 
vanish and the corresponding dielectric constant will become infinite. 
The figure shows that there is a transition at about 1400 atmospheres 
and the trend of the curve suggests that there will be a further transition 
at sufficiently high pressures. Bancroft made similar sets of measurements 
at other temperatures, but his results are best shown by drawing contours 
of constant y,,;* on a pressure-temperature diagram. This is done in 
fig. 2, and on the same diagram are drawn lines of constant volume. 
This diagram shows clearly that two transitions take place under 
conditions of constant volume as well as constant pressure, and that the 
temperature range of the spontaneously polarized state increases as the 
pressure increases. The values of reciprocal susceptibility plotted in 
fig. 2 all refer to a state of zero polarization. In the spontaneously 
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Contours of constant reciprocal susceptibility on a pressure-temperature 
diagram (Bancroft 1938). 


polarized state, however, Bancroft’s measurements refer, of course, to 
crystals with the stable spontaneous polarization. The values for zero 
polarization can be deduced as follows. At a given temperature and 
pressure the free energy can be expanded in even powers of P according 
to the equation 

G=G,.+4xiPY+BPy, ie ie a 22.) 
where higher powers of P, can be neglected if P, is small. The electric 
field H is given by 


0G : A 23 
i= OP, === GT P,+4BP, > fpr) ee ee (2. ) 

and the reciprocal susceptibility by 
ea 2250 SCTE a re C2 


oP, 
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A spontaneously polarized state occurs when X11* is negative. Then if 
the field is zero we have 

O= 177° P,+4BP{, 
and hence P, is zero or 

PP=—x11°/4B. 
It can easily be shown that the last solution is the stable one and if we 
substitute in (2.4), we have 

0L 

aP, ia 2x ie 
that is 


ak, 


X11 =— tap (2.5) 


Thus we have the value of y,,* in terms of the reciprocal susceptibility 
for the spontaneously polarized state. This equation was used in plotting 
the contours in fig. 2. 

When we consider the properties of an atomic model it is convenient 
to take the polarization and linear dimensions of the model as variables, 
so that the most directly derived quantity is the reciprocal susceptibility 
at constant strain y,,”. This is related to y,,* by the formula 
(Devonshire 1954, eqns. (2.4)) 


Xa =X +O Cas" 2p) eee (ZG) 
In this equation b,, is the piezoelectric constant that gives the ratio of 
shear strain in the yz plane to polarization in the z-direction, and c¢,,” 
is the ordinary shear elastic constant measured at constant polarization. 
The quantities b,, and c,,’ have no anomalies at the transitions and 
do not vary much over the temperature range we are interested in; 
they can be determined by piezoelectric measurements (Mueller 1940). 
Alternatively they can be determined directly by making measurements 
at such high frequencies that the piezoelectric response does not take 
place and then compared with y,,* (Mason 1947). In either case we 
find that 
X17 — X11 0-04, Oh 


and only varies slightly with temperature (Depee 1954, p. 97). 
It is not known how this quantity varies with pressure, but probably 
this variation is not large. Since the vanishing of y,,” is the condition 
for a second-order transition under conditions of constant strain instead 
of the normal conditions of constant stress, and y,,” is always larger than 
X11°, it follows that if the crystal could be clamped so as to prevent 
spontaneous strain then the range of the spontaneously polarized state 
would be reduced. In fact the range of the spontaneously polarized 
state of the ‘clamped’ crystal on the pressure-temperature diagram 
will be that within the contour y,,*=—0-04. This means that at zero 
pressure the temperature range of the spontaneously polarized state will 
be just about reduced to zero, though at higher pressures there will be 
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still a finite range. This was first pointed out by Mueller (1940). It 
does not seem possible to test this prediction directly ; for although the 
erystal can be clamped in such a way as to prevent the strain that 
accompanies induced polarization this will not affect the spontaneous 
polarization ; the reason is that the spontaneous polarization normally 
takes place in domains arranged in such a way that the net strain over 
the whole crystal is small. 


§ 3. Mason’s THEORY 
Mason’s theory of Rochelle salt was based on the crystal structure as 
determined by Beevers and Hughes (1941). The x-rays did not show the 
positions of the hydrogen atoms, but the positions of the oxygen atoms 
appeared to indicate that some of the hydrogen atoms could migrate 
between two positions of equilibrium as shown in fig. 3; the figure is 


Fig. 3 


asymmetrical because the atoms have different binding energy in the 
two positions. A movement of a hydrogen atom from one position to 
another will change the electric moment of the cell, and the system may 
be considered as a dipole that can reverse its sign, and has a different 
energy for the positive and negative orientations. The movable hydrogen 
atoms occur in pairs in the unit cell with opposite orientations, so that 
when both are in their positions of lowest energy there is no net dipole 
moment. So far we have only considered the isolated unit cell. When 
the substance is polarized, however, there will be an electrical interaction 
between the dipoles in different cells, which will lower the energy. At 
very low temperatures the atoms will arrange themselves so as to make 
the total energy a minimum. This arrangement may be either the 
‘antipolar ’ state in which each hydrogen atom is in its position of 
lowest energy, so that half the dipoles point each way and there is no 
polarization—or the polar state in which all the dipoles point the same 
way ; it depends on whether the lowering of energy due to the interaction 
of the dipoles is sufficient to compensate for the increase of energy due 
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to half the hydrogen atoms being in the higher energy state. At high 
temperatures the atoms will be more or less randomly arranged, as this 
state has the largest entropy. If conditions were such that the substance 
passed successively through the antipolar, polar and non-polar random 
states, then the existence of two transition temperatures would be 
satisfactorily accounted for. 

Recent work on the structure of Rochelle salt by neutron-diffraction 
by Frazer etal. (1954) has shown that the structure determined by 
Beevers and Hughes was not quite correct. The particular hydrogen 
atoms that were supposed to move clearly could not do so. They showed, 
however, that certain others had two possible positions of equilibrium. 
Thus the essential features of Mason’s model are not altered. His 
calculations of energy barriers, dipole moments, etc., have no significance 
now; however, in this paper we are not concerned with these details, 
but only with the general features of the model, and whether it can give 
the double transition observed in Rochelle salt. The model may be 
applicable to other substances. 

The movement of the hydrogen atoms is not, of course, the only 
source of polarizability in the material. There is also the electronic 
polarizability, but we shall see later that this does not affect the main 
conclusions, so for the moment we shall not take it into account. 

Let there be NV dipoles per unit volume of each group of fixed moments 
+p, and let the difference in energy between the two dipole positions 
be 24. At any time let 3N(1-+<2,) of the first group and }N(1+2,) of 
the second group point in the positive direction. Then the polarization P 
will be given by 

P=Nupu(x,+2,), Sn ee een oe 
and the dipole interaction energy will be —}yP?, where y is a constant 
depending on the configuration of the dipoles. If the dipoles of the first 
group have their minimum energy when pointing in the positive direction 
then the energy of the system per unit volume (referred to the antipolar 
state as zero) is given by 

U=—FyN*u?(e1+%_)?+NA(2+e,—2,). . . . (8.2) 
The entropy of the system is given by the usual formula in terms of the 
number of configurations 
ee N! ; N! 
SSE NE IN ta) Naa)! + * 8p iN(—a,) | 

= —kN{(1+2,) log (14-a,)-+(1—a,) log (1—a,)+(1+a,) log (1-++2,) 

+ (129) log (1—a_)—4 log 2}, . . (3.8) 
and since the Helmholtz free energy A is U—T'S this is given by 
A|NkT = —}a(%4+%2)?+B(2+2,—a,) 

+3{(1+2,) log (1+-a,)+(] =n 


log (1—a,)-+(1++a,) log (1+2,) 
log (l—2,)—4 log 2}, . (3.4) 


where a=yNwkT, B= A/kT, a eed 
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At equilibrium A is a minimum with respect to both x, and x, and hence 


1 
0=—a(x,+a)—B+3 log ae > 
1 
1+ 2, 


=—a(21+2%.)+8+3 log ; 
so that 
2,=tanh {a(7,+2,)-+ f}, 3.6 
vo—=tanh {a(%,+2,)—f}, leo, 
and therefore 
+ ,=tanh {a(7,+a,)+}+tanh {a«(x,+a,)—P 
— 2 sinh 2a (%,+2,) 
~ cosh 2a(%,+2:5)-+-cosh 28 * 

If we plot the right-hand side of eqn. (3.7) as a function of x,+2, 
it takes two forms according to whether cosh 28 is less than or greater 
than 2, as shown in fig. 4. In the first case there will be a solution other 
than zero if the slope at the origin is greater than 1. This will be so if 


(3.7) 


4a 
ee I-F cosh 2B" (3.8) 
Fig. 4 
2.0 
1.0 


Xi+ x,— 


The transition temperature will be given by 

4~%—=1-++cosh 28, nes got eee (3.9) 
and the transition will be a second order one. In the second case the 
transition will be a first order one and will take place when a line through 
the origin with unit slope cuts off equal areas. 

We have now to consider whether eqn. (3.7) can have two solutions 
for different temperatures. If the volume is kept constant as the 
temperature changes it seems reasonable to assume that the quantities 4 
and yNp2 do not change. We can therefore write 

B= 0a, rere. (SO) 


where 0= A/yNp?, ss ek (StL) 
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and so is independent of temperature. Now at the absolute zero the 
substance will be either in the completely ordered polar state (x, =x,=1), 
or in the completely ordered antipolar state (vj;=—a,=1). In the first 
case the free energy is —2«-+28; in the latter case it is 0. Thus if the 
antipolar state is to be the stable one at low temperatures we must 
have @>+1. If this is go, however, can we easily show that eqn. (3.9) 
has no solution, so that no second-order transition to a polar state is 
possible. It can also be shown that no first-order transition is possible. 

Thus Mason’s model as it stands does not give two transition temperatures 
at constant volume. We shall now show, however, that with a slight 
modification it will do so. An obvious modification to the simple model 
we have been using is to assume that the two states of the hydrogen atom 
have different weights. Classically this corresponds to the fact that the 
atom may have different vibration frequencies in its two equilibrium 
positions. One would expect it to have a lower frequency in the position 
with the higher energy, so that this state should have the greater weight. 

Let us suppose that the higher energy state has a greater weight by a 
factor aw (it is clearly only the ratio of the weights of the two states that 
matters). Then the energy of the system will be unaltered, but the 
entropy will be increased by an amount $Nk(2+%7,—x,) log w. Equation 
(3.4) for the free energy retains the same form, but instead of (3.5) we 
have _ yp? 

Paley 


= ia eM MEER 9151055) 
where, as before 0= A/yNp?, 
@=¢ log! 47s 15 ieee 5 oe er) 
Equation (3.9) now becomes 
4a—=1-+cosh 2(@a—¢). x 8 STS) 


Unless ¢ is very small this equation in « will, in general, have two solutions 
for suitable values of 6. In fig. 5 we show @ plotted against 1/x. Since 
I/« is proportional to 7 this figure shows that in a suitable range of 6 
we have two transition temperatures, and that the upper one is always 
a second-order transition, but the lower one may be first or second order 
according to the value of 0. 

We must now consider the effect of the electronic polarization, which 
we have neglected so far. If we denote this by P, there will be two 
additional terms in the energy and hence in the free energy, namely 


pOP 2—bP Nu(x,-Ea.): Vas, sie See (Oho) 


The second term represents the interaction energy between the electronic 
polarization and the dipole polarization, The first term is the net energy 
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required to produce the electronic polarization; the constant in the 
term can be determined in terms of the refractive index, which depends 
on the electron polarizability alone. Since A has to be a minimum with 
respect to P, as well as to x, and x, we have an additional equation, 
namely 

0=aP,—bNp(x,+2,), eS he i ee (PLO) 


and if we substitute for P, in the expression (3.15) then the additional 
term in the free energy becomes 


—3(b7/a)N*u?(x4-+-25)?. 
Thus comparing it with (3.2) we see that the sole effect of the electronic 
polarization is to replace y by y’ where y’=y-+(b?/a). 


Fig. 5 


1-8 1.6 [.4 (iS 
Curves at which a transition takes place. The full lines indicate a second order 
transition and the dotted lines a first order one. 


Let us now consider the dielectric susceptibility. We have 
0A cannes 04 


oP Wa oe Ore 
and hence 
E=aP ,—bNp(«4,+2,) er Oe chil hte nl ache (ost as) 
” l is 
=—bP 4 “| a(v,+x,)—B+4 log a 2 (orks 0) 
als ey 
a 1 
=P 4 Tf ale ta) tittle}. 8.179) 
‘ates a" s 


from (3.4) and (3.15). The total polarization P is given by 
P=P ,+Np(x,+%,). fee ee 015-18) 
From (3.176) and (3.17¢) we have 


2,=tanh lz (HBP, + alary+t9) +8. 


%_=tanh { Fo-+0P,)+aletre)—B} 
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¢,+2,—=tanh a W+bP baler tan) +P} 


+tanh la HOP.) +a(e,+02)—B| 
2 sin 2f{u/kT(E+bP,)+a(4%1+%2)} 
~ cosh 2{u/kT (H-ObP,)+a(a,+a5)}--cosh 28° 
From (3.17), (3.18) and (3.19) we can determine # in terms of P, and 
hence find the dielectric constant. If we substitute from (3.17a@) and 
(3.18) in (3.19) we find that 


2 2 okeT Pak 
Me 2 sinh 7 enon {(a+2b)E+b?P+ NE (a )} 
(atbd)Nu — Qu : akT } 
———— —— (aP—H sh 2 
cosh i (ab) {(a-+ 2b) H-+-b2P + Me (a )}+ cosh 28 
If # and P are both small then we have 
4Nu2 
(1+ cosh 28)(aP—E)= iF { (a+ 2b)E+6?} +-40(aP—E), 

and hence 


EH _ a{i+cosh 2B—40}—(4Ny2/kT’) b? 


= : . =. (3.20 

P 1+ cosh 28—4a+ (4N2/kT) (a+ 26) (3:20) 
This can be further simplified if we are fairly near the transition 
temperatures. 


Here H/P will be small and we shall have approximately 
a{1+cosh 28—4a} ~4Np?/b?2kT. 
This can be substituted in the denominator, and then we have 
1 k 
= Somer 


Et 
pede ts Se Re Ay = 
(Lb aN {1+ cosh 28—4« } 


ab? a 
aqpp + 822) 
If we substitute for « and £ from (3.5) then we have 
EK ane kT a*y+ab? : 
P= GEbp ane (i teosh 28)— Cena wo ER) 
where 
B= 3 logo, 
that is 
12S all shar gna pyplceae 
P Gpabare L osh 28 }—y | 
a! 
Gia [1+ cosh 28—4«’], ) herleie (ores) 
where a =y’ NurVkT. 


(3.24) 
§ 4. COMPARISON With EXPERIMENT 


We must now compare the predictions of the theory with the experi- 
mental results of Bancroft, 


The simplest comparison is that of the 


(3.19) 
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transition temperatures. The theoretical calculations refer to transitions 
at constant volume ; but the experiments were made at constant pressure. 
The transitions at constant volume (if they are second order) occur when 
X11" 1s zero and hence, from (2.7), when y,,* is —0-04. Thus the experi- 
mental contour of y,,* equal to —0-04 plotted in fig. 2 should correspond 
to the theoretical contour of the transition. Now the theoretical contour 
of the transition depends on three parameters y’Np2/k, 4/k, 6, which we 
may reasonably assume to be functions of volume only. There are two 
transition temperatures for each volume; thus if we knew one of the 
parameters as a function of volume we could determine the other two 
by a comparison of theoretical and experimental results. We find that 
the values of the other parameters are not sensitive to the value of ¢, 
and this is unlikely to vary much with volume. We shall therefore 
arbitrarily take ¢ to be 0-6. We then find the following values for the 
other parameters as functions of volume (table 1). If we take a different 
value of ¢, say 0-3, then we get the following set of values (table 2). 


Table 1 


volume (cm3/mole) | y’Ny?/k Alk 


156-2 180° oll 
1543 193° 332° 


Table 2 


volume (cm3/mole) | y’Nu2/k | A/k 


156-2 203° 260° 
154-3 218° Qh be 


Thus in each case we get a similar variation with volume. Over the 
small range of variation considered it may be taken to be linear. We see 
that both parameters increase as the volume decreases. Now there is 
no particular a priori reason for predicting how 4 should vary with the 
volume, but it seems surprising that y’Ny? should increase as the volume 
decreases. Of the factors involved N must, of course, increase, but the 
change is very small; on the other hand ,. must surely decrease as the 
two equilibrium positions of the hydrogen atom come closer together, 
and this will far outweigh the change in N. The increase in y’Nu® must 
presumably therefore be due to an increase in y’. Now we have 


y=y+h/a. ea, ee (ART) 


Of the two terms in y’, y depends only on the arrangement of the dipoles, 
and there is no particular reason why it should change with volume. 
The same applies to b, which determines the interaction of the dipoles 


4A2 
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and electronic polarization. The value of a determines the electronic 
polarizability and can, in fact, be calculated from the refractive index ; 
we should expect the electronic polarizability to increase somewhat as 
the volume is diminished, but not so much as to cause the observed 
change in the parameter. In fact we should expect the main change 
with volume to be in the parameter dipole moment , and this would be 
in the wrong direction. 

We can also gain information concerning the constants in the equation 
by comparing the predicted polarization and dielectric constant with 
the experimental values. Now the dipole polarization is Nu(x,+7,), and 
to this must be added the electronic polarization P,. The relation 
between them is given by (3.16). Thus the total polarization is given by 


P=Np(«,+x,)(1+0/a). no Bo are eee) 


Now the maximum value of the polarization at zero pressure is 740 e.s.u., 
and the corresponding calculated value of x,-++-x, is 0-4. Hence we have 


Nu {1+ (b/a)}=1850, 2. ones Setter 
and JN is 3:81 x 1021, so that 
u(1-+b/a)=0-49 x 10-18 e.s.u. eg oe Ste 


The value of the dielectric susceptibility is given by eqn. (3.23). Now 
above and below the transition temperatures the reciprocal susceptibility 
is nearly a linear function of temperature, so the simplest quantity to 
compare is the slope of this. According to Bancroft we have at the 
upper critical temperature 

0 (E* 


—) =65: —8 
at\P), 50x 1079, 
a =} =—0-061x10- kg/cm? 
op\P) oe x per Kg/cm*. 


The theoretical values will refer to rate of change at constant volume, 
but this can be deduced from the above values if we use the values given 
by Bancroft for thermal expansion and compressibility. From these 
we can deduce that 
a 
(sr). = 23-8 kg/cm2/deg, 
and hence that 


g E —(5-6 1°45° =3 5 es 
Ale i —1-45)10-8=4-15x 10-3, =. . . (4.5) 
and 
0 (E 
i —]-9 
Pa(p) =1 23. . a rere 
We can calculate that at the upper critical temperature 


r Cael ; S ’ 
Ya a + cosh 2B—4« )} =o, 


Theory of Rochelle Salt 1039 


and hence from (3.23) we have 


a*y’ 
or 
a ese (bia et 2 he Oe eer) 


By combining the values given in table 1, however, with the values 
given in (4.3) and (4.4) we can deduce that 

Pere Te Vee (ee Verges df Niied Guaeve 1(4E8) 
Thus the values of our constants calculated from different experimental 
results are not entirely consistent. 

The quantity @ can be deduced from the refractive index, since it gives 
the reciprocal value of the electronic susceptibility, and this determines 
the refractive index. The refractive index is 1-49, and hence the corre- 
sponding dielectric constant is (1:49)?, that is 2-22, and this is related 
to a by the formula 

2-22—14+4z/a, 

so that 

a=10-3. 
The value of 6 can only be guessed ; it is the factor giving the interaction 
of two sets of dipoles, and if these were at random its value would be 
47/3. These values of b and a would make y’ equal to 16 according to (4.7) 
or to 54 according to (4.8). The values of y given by (4.1) would be only 
a little less. But y is of the same nature as b, since it determines the 
interaction between the dipoles formed by shift of the hydrogen ions, 
and hence these values are rather high. 

Thus we see that Mason’s theory, when modified, can account fairly 
well for the observed facts, but there are still some difficulties. Some of 
these could be removed by further refinements. It seems difficult, 
however, to account in a reasonable way for the fact that the range of 
transition temperature increases with decreasing volume. As we have 
just seen this seems to imply an effective dipole moment that increases 
with decreasing volume, and this is difficult to explain. 
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ABSTRACT 


On January 24, 1956 observations were made at a wavelength of 
3-7 metres, of the lunar occultation of the Crab Nebula. A comparison 
of the time of obscuration of the nebula with the computed value enabled 
the refraction occurring in the lunar atmosphere to be estimated. An 
electron density at the moon’s surface of about 10° cm~* was derived 
for cases of both an atmosphere in hydrostatic equilibrium and for an 
atmosphere which is continuously escaping ; a figure which corresponds 
to a surface density of the lunar atmosphere about 2 10-1 of that of 
the earth’s atmosphere at normal temperature and pressure. 


§ 1. INTRODUCTION 

OpTicaL measurements have failed to detect an atmosphere surrounding 
the moon; the most sensitive observation was that of Dollfus (1952), 
who, by measuring the polarization of the scattered light, showed that 
the density of the lunar atmosphere is less than 10~° of the density of 
the terrestrial atmosphere at sea level. A density of 10~-® of the terrestrial 
atmosphere derived by Lispkij (1953), has been proved to be unacceptable 
by Opik (1955). 

It has been shown by Link (1956) that the density of the lunar 
atmosphere may be derived from measurements of refraction of radio 
waves at the sunlit limb of the moon, when a radio star is occulted. 
Preliminary results obtained from observations of the occultation of the 
large diameter radio source in the constellation of Gemini (Elsmore and 
Whitfield 1955), have indicated that the density is less than 10-12 of the 
density of the terrestrial atmosphere. This paper is concerned with 
deductions from the occultation of another radio source, the Crab Nebula, 
which was observed at Cambridge on January 24, 1956 (Costain, et al. 
1956). In § 2 the radio observations are described and their interpretation, 
in terms of the electron density surrounding the moon, is discussed in 
§3. In §4 an attempt is made to deduce the density of the lunar 
atmosphere from the measured electron density. 


{ Communicated by the Author. 
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§ 2. THe Rapio OBSERVATIONS 


The occultation of the Crab Nebula was observed at wavelengths of 7-9 
and 3-7 metres. The measurements at 7-9m were seriously affected 
by man-made interference and it was only possible to observe the 
immersion, when the source disappeared behind the moon. At a 
wavelength of 3-7m accurate observations were made throughout the 
occultation, from which a curve was plotted of the observed intensity 
of radio noise from the nebula. As the diameter of the source at this 
wavelength is approximately 5’ of arc, the intensity, at immersion, 
decreased to zero in about 10 minutes, and was restored about an hour 
later at emersion, when the source reappeared at the opposite limb of 
the moon. Since the moon’s centre passed near the centre of the source, 
different points on the nebula were covered and uncovered in the same 
sequence. The occultation curves obtained during emersion and 
immersion were therefore complementary and corresponding points 
on the curves were found to be separated by an interval of 59-6+-0-26 
minutes. The error represents the scatter of times derived from different 
parts of the occultation curve. 


§ 3. INTERPRETATION 


The time of obscuration, in the absence of refraction, was calculated 
from elements of the occultation predicted using Bessel’s method and was 
checked by interpolating topocentric coordinates, kindly provided by 
Dr. Sadler of the Nautical Almanac Office. The calculated time, 
59-21 minutes, is 0-4--0-26 minutes less than the observed time. 

Although the difference between the observed and calculated times 
represents a barely significant result, it will be assumed for the present 
to have a meaning in deriving a value for the electron density ; the 
observations, in any case, provide a figure for the maximum permissible 
density, as will be shown later. 

Two phenomena may contribute to this time difference ; diffraction 
at the limbs of the moon (§ 3.1), and refraction in the lunar atmosphere 
(§ 3.2). 

3.1. Diffraction 

The departure from a straight edge, due to both the curvature of the 
moon’s limb and the mountainous surface is small compared with the 
size of a Fresnel zone. It was further found, that for the present case 
of an extended source, the effect of Fresnel diffraction is negligible. 


3.2. Refraction 

On the sunlit side of the moon, it is likely that the atmosphere, which 

had been exposed to the sun’s radiation continuously for Il days, was 
ionized. The refractive index, », of an ionized gas is given by 

a ve Sisk 10% 
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where N is the electron density and f is the frequency in Me/s. Since 
p==1, we can write 
4-05 x 10-6 


1p 


The angle of refraction « of the radio waves will depend upon N and 
its gradient and two different types of distribution will be considered. 
Firstly, an exponential variation of electron density, which would be 
established if the atmosphere were in hydrostatic equilibrium, and 
secondly, an inverse power law, which would occur if the atmosphere 
were continuously escaping from the lunar surface. 


3.2.1. Exponential variation of electron density 


On the assumption that the Chapman formula is applicable, Link (1956) 
has shown that the electron density at a height h above the moon’s 
surface is given by 

N=N, exp (—A/2H) 
where N, is the electron density at the surface, and H the scale height. 
For this distribution, the angle of refraction «, which is assumed to be 
small, has been calculated by Link to be 


a= (1—p)/(wa/H) 


where a@ is the moon’s radius. Substituting for 1—y from above 


DA06x 102/78) h 
r= Tiare 7) =%0 &XP | — a9 


where a, is the angle of refraction for a ray grazing the moon’s surface. 


3.2.2. Power law variation of electron density 
If the electron density at a height h is given by 


where / is measured in units of the moon’s radius, and y is a constant, 
the angle of refraction is 

es 4:05 10-5 2Noy f° 

jet EUs 

where 7 is the angle of incidence of a ray entering a spherically stratified 
medium. If the mean vertical component of velocity of the ionized 
gases is equal to the velocity of escape, then y=3/2; but if the velocity 
is appreciably greater than the escape velocity, then y=2. 

The angle of refraction of rays which graze the surface of the moon 
is given in the table for different values and distributions of electron 
density. Since the moon moved its own diameter, which subtended 
32’ arc at the earth, in 59 minutes, the observed increase in the time of 
obscuration corresponds to «,—13-4"-+-8-7" of arc, if it is assumed that 
refraction occurred only at the sunlit limb, and half this value if it 
occurred equally at both limbs. It can be seen that the electron density 


sin’ ¢ da 
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lies in the range 10? to 104 cm-® for y=3/2, y=2 or for scale heights in 
the region 5000 km to 80km. If the scale height lies in the region 
80 km to 20 km, the electron density is about 10? cm-2; under these 
circumstances the extent of the refracting region necessary to produce 
the observed increase in the time of obscuration determines the lower 


limit of the scale height. 


Table giving the Angles of Refraction, %) for a Grazing Ray 


Electron density at surface 
(=, electrons cm?) 


10° 10+ 10 


H=5000 km Zz 2 13-2” 1-3" 
H= 500 km 6’ 55” 41-5” 4-2" 
H= 100 km 15’ 30” 33" TDS se 
H= 50km 34’ 40" 3’ 28" 20:8” 
y=3/2 4’ 23” 26:3” 2-6” 
y=2 6’ 36” 39-6" 4-0" 


§ 4. Discussion 


The observations provide a measure of the electron density near the 
moon’s surface, in excess of that of the surrounding medium. Before 
relating these data to the two types of atmosphere it is necessary to 
consider the possibility of the density being increased on one side of the 
moon, due to its motion through the interplanetary matter. Observations 
of the zodiacal light (Siedentopf et al. 1953) and of ‘ Whistler ’ atmos- 
pherics (Storey 1953) have indicated that the electron density near the 
earth is approximately 600 cm~*, hence the orbital velocity of the moon 
round the sun (30 km sec~!) might be sufficient to cause an appreciable 
increase in density near the leading face of the moon. 

The magnitude of the increase will depend on the surface effects, 
when the proton-electron gas, having a temperature of the order 
of 10000°K, comes into contact with the lunar surface at a temperature 
of 160°K. It seems probable that in such an encounter, the majority 
of the particles will recombine at the surface and will subsequently 
escape before being reionized. Since the observed effects can only be 
explained if all the particles, after collision, leave the surface with a 
velocity of 2-0 km sec~!, corresponding to the temperature of the dark 
surface, and retain their charge, it seems unlikely that this mechanism 
is important. 

The electron density is therefore probably caused by the action of 
solar radiation on the lunar atmosphere and estimates of the gas density 
will now be made for the cases of a ‘permanent’ and an escaping 
atmosphere. 
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4.1. Atmosphere in Hydrostatic Equilibrium 


It is difficult to estimate the degree of ionization of the lunar atmosphere 
by comparison with the terrestrial ionosphere, because its composition 
probably differs from that of the earth’s upper atmosphere in the 
following ways. 


4.1.1. Chemical composition 

Whereas the earth’s atmosphere above 400 km consists almost entirely 
of N, and O, the atoms and molecules comprising the ‘ permanent ’ 
lunar atmosphere must have a molecular weight greater than 42 if the 
temperature is equal to, or greater than the surface temperature of 
380°K ; otherwise the molecular velocities will exceed the velocity of 
escape. 


4.1.2. Temperature 

Both the mechanism responsible for heating the earth’s upper atmos- 
phere, and the temperature above a height of 200 km are uncertain. 
By extrapolating results obtained from observations using rockets, Bates 
(1954) has derived, for two different models, the values of 900°K and 
1700°K for the temperature at 500 km. Nicolet (1956) and Urey (1957) 
have shown that even higher temperatures must be assumed to account 
for the observed rates of escape of ?3He and 4He. However, the 
temperature of the lower region of the lunar atmosphere is probably 
not very different from that of the surface (380°K on the sunlit side), 
and it is difficult to estimate how the temperature varies with height, 
when the mechanism for heating the terrestrial ionosphere is, as yet, 
not understood. 

If the temperature is 400°K, the scale height of the ‘ permanent ’ 
lunar atmosphere would be about 50 km and hence, it can be seen from 
the table that the electron density required to produce the observed 
angle of refraction is 103 em-. 

In the terrestrial ionosphere above a height of 400 km, about 1/100 
of the total number of particles are ionized. Ags the molecules in the 
‘permanent ’ lunar atmosphere have a high molecular weight, they are 
more readily ionized than those in the ionosphere. However, diffusion 
to the surface of the moon, where the charged particles will recombine, 
may tend to reduce the proportion of charged particles. It is difficult 
to make a quantitative estimate of the latter effect, but provided it is 
not of undue importance, it is probable that at least one in a thousand 
particles is ionized. The total number of particles in the lunar atmosphere 
is therefore unlikely to exceed 10°cm-*; the density of the lunar 
atmosphere would then be 10-® of that of the earth’s atmosphere at 
sea level. 

If it is considered that the observations provide only an upper limit 
to the electron density surrounding the moon, then the density of the 
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‘permanent’ lunar atmosphere must be less than 5x 10-13 of that of 
the earth’s atmosphere at sea level. 

It is unlikely that the atmosphere will condense on the surface at the 
dark side of the moon, for, although the surface temperature is only 
about 160°K, the lunar atmospheric pressure is very much less than that 
required to liquify the probable constituents. 


4.2. Hscaping Atmosphere 


It is probable that certain gases, of such low densities that hitherto 
they have been neglected, are continuously escaping and are being 
replaced at the moon’s surface. Such mechanisms responsible for this 
are (1) radioactive decay of part of the rocks comprising the moon’s 
crust, leading to the production of helium and argon for example, 
(2) vaporization of a fraction of the 210° g of meteors which strike 
the surface of the moon daily, and (3) release of gases from the interior 
of the moon. 

Provided that about one in a thousand of the escaping molecules is 
ionized, the density of the atmosphere at the surface of the moon is 
about 5x 107-18 of the earth’s atmosphere at normal temperature and 
pressure. 


§ 5. CONCLUSION 


The observations provide no indication of the relative abundance of 
the two forms of atmosphere ; such information could only be derived 
from detailed measurements of the variation of refraction with distance 
from the moon’s limb. In the event of approximately equal abundances, 
the density of the atmosphere at the surface of the moon is about 
2x 10-18 of that of the earth’s atmosphere at sea level, and is unlikely 
to exceed 6x10- of that of the terrestrial atmosphere at normal 
temperature and pressure. 

As the mean free path in the lunar atmosphere is very long, it is probable 
that the ionization may extend over the dark side of the moon. In the 
extreme case of equal refraction occurring at both limbs of the moon, 
the above values of the density are divided by two. 
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§ 1. INTRODUCTION 


Two techniques have recently been developed which allow dislocations in a 
number of transparent inorganic crystals to be studied with the micro- 
scope. The first is based upon the fact that etching solutions often attack 
particular crystal surfaces more rapidly at points where dislocations emerge 
than elsewhere. Much valuble information about the properties and be- 
haviour of the dislocations may then be derived from the study of the etch 
pits. The second depends upon the internal separation of a decorating 
substance along the dislocation lines. 

The most useful observations will be made when it is possible to use the 
two methods in succession for the study of the same dislocations ; this has 
so far been possible only in the case of silicon. Here, Dash (1956, 1957) 
has appled both techniques to the same crystals and established the 
correspondence between etch pits and decorated dislocations. It has not 
yet proved possible to decorate the dislocations within crystals of lithium 
fluoride but Gilman and Johnston (1956, 1957) have made many observ- 
ations of the motion of dislocations in the crystals with the etch pit 
techniques which they have developed. Barber et al. (1957) have worked 
out a useful method for decorating the dislocations in crystals of sodium 
and potassium halides with gold but reliable etching of the dislocations in 
the crystals prior to decoration has not yet been achieved. In all these 
inorganic crystals, the properties of the dislocations which are always 
present in the annealed crystals are different from those of the additional 
dislocations which appear during deformation. 

The purpose of this note is to describe the etching of both types of dis- 
locations in annealed and plastically deformed crystals of silver halides. 
The correspondence between the etch pits and the points of emergence of 
dislocations has been established by subsequent exposure which causes 
silver to separate along the dislocation lines (Hedges and Mitchell 1953). 


§ 2. EXPERIMENTAL METHODS 
The silver chloride, silver chloro-bromide, and silver bromide sheets 
which are being used for this work are grown from the melt by the method 
described by Clark and Mitchell (1956). After separation from the glass 
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plates, they are annealed for 8 hours in an atmosphere of the halogen at a 
temperature between 200 and 400°c. The halogen is then condensed by 
immersing a side arm in liquid air for | hour before the furnace is allowed 
to cool. The crystals are deformed by elongation and by indentation with 
hemispherical glass indenters and then etched and exposed to reveal the 
dislocations. . 

Although many solvents attack the surfaces of silver halide crystals 
preferentially around the points of emergence of dislocation lines, ex- 
perience has shown that a 3N solution of sodium thiosulphate provides the 
most reliable differential etchant. The solution is allowed to act for the 
shortest possible time ; the specimens are simply immersed in the solution 
and transferred immediately to water. The etch pits may subsequently 
be enlarged, usually with some loss of resolution, by treatment with a 
much weaker solution of sodium thiosulphate which does by itself not 
etch the dislocations reliably. 

After they have been thoroughly washed and dried, the etched crystals 
are exposed to the unfiltered radiation from a 250 watt high pressure 
mercury vapour lamp for 30 seconds at a distance of about 30 cm from the 
lamp; this causes the separation of silver along the dislocation lines to a 
depth of at least 40 microns below the surface. 


§ 3. EXPERIMENTAL OBSERVATIONS 


The thin sheets which are produced by crystallizing discs of molten 
silver halides between glass plates usually consist of a few large crystals with 
surfaces almost parallel to (001) planes. After annealing, the presence 
of dislocations in crystals with precise (001) surfaces, whether the crystals 
have been plastically deformed or not, cannot usually be established either 
by etching or by the separation of internal silver (fig. 1, Pl. 39). These 
crystals do not show sub-structures and have a very low reactivity towards 
photographic sensitizers and developers (Evans et al. 1955). Very few 
dislocations are observed in crystals with surfaces almost parallel to (001) 
planes (fig. 2, Pl. 39). As the deviation increases, a few isolated walls first: 
appear (fig. 3, Pl. 40) followed by the well defined sub-structures which have 
been studied in detail in our previous work (Mitchell 1957). The disloca- 
tions thus etch and due to subsequent exposure to light are decorated by 
the internal separation of silver to a steadily increasing extent as the angle 
between the surface and the (001) planes increases from 0° to about 10°, 
the maximum deviation usually found in these crystals. 

Tt will be observed in the photomicrographs that the etch pits which are 
produced by the 3N solution of sodium thiosulphate have smoothly curved 
sides and rounded bases. They may be enlarged to pyramidal pits with 
apparent crystallographic facets by subsequent immersion in a weaker 
solution of sodium thiosulphate. 

A careful study of the correspondence between etch pits and disloca- 
tions has shown that not all etch pits are associated with dislocations 
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which may be made visible by the separation of internal silver. Neither 
do etch pits appear at all the points where dislocations, which can be made 
visible during exposure, meet the surface. The correspondence depends on 
the precise crystallographic orientation of the surface and upon particular 
properties of individual crystals; it is complete in many crystals but 
absent in others. 

The dislocations which appear during plastic deformation at room 
temperature give smaller etch pits and print out to a smaller extent than 
those which are already present in annealed specimens. The two types of 
etch pit are shown in fig. 4, Pl. 40. This specimen was etched, deformed with 
a hemispherical indenter, and then re-etched. The larger pits which are 
not as well defined as in figs. 2 and 3 are associated with the grown-in 
dislocations and the smaller sharply defined pits with dislocations which 
appeared during plastic deformation. The most sensitive conditions for 
the study of the behaviour of dislocations during the deformation of silver 
chloride crystals occur when the crystals have surfaces almost parallel to 
(001) planes ; in these crystals, the apparently small numbers of grown-in 
dislocations can only just be rendered visible by etching. The small and 
shallow etch pits corresponding with dislocations on slip bands produced 
by the deformation of such crystals can easily be resolved with the phase 
contrast microscope. 

With exposure to light silver separates along the dislocation lines which 
appear during plastic deformation to form an internal latent image, and 
this accounts for the greatly enhanced internal sensitivity of such crystals. 
The internal image may be developed after the removal of the surface 
with a dilute solution of potassium cyanide. In favourable circum- 
stances, the exposure may be continued until the dislocation lines become 
visible within the crystals (Mitchell 1957). 

Our observations indicate that the dislocations in the annealed crystals 
which give rise to deep etch pits and along which visible amounts of silver 
separate rapidly during exposure do not play an important role during the 
early stages of plastic deformation. It may be concluded that their 
mobility is reduced by the segregation of impurities or of vacancies during 
cooling to room temperature after crystallization or during subsequent 
annealing. When the crystals are deformed, dislocations appear within 
the elements of the substructures as in the experiments on lithium fluoride 
of Gilman and Johnston (1956, 1957). 


§ 4. CONCLUSION 


The silver halides provide the only inorganic crystals, which have 80 far 
been studied, in which the presence and properties of the dislocations 
associated with both annealing and deformation structures can be estab- 
lished by the application of etching and decoration: techniques at room 
temperature. The dislocations can be rendered visible under the actual 
conditions of deformation and modifications due to the decoration 
techniques are therefore excluded, 
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A Demonstration of Parity Non-Conservation in B-Decay 


By M. A. Gracz, C. E. Jonnson, R. G. Scurntock and C. V. Sowrer 
Clarendon Laboratory, Oxford 


[Received June 17, 1957] 


Since the suggestion of Lee and Yang (1956) that parity conservation 
might not hold in weak interactions, attention in recent work on nuclear 
orientation has been focused on the angular distribution of B-particles 
emitted from polarized nuclei. We have established the feasibility 
of detecting £-particles by the blackening of a photographic film inside a 
cryostat at a temperature of 1°k and have used this method of recording 
to measure the angular distribution of the B-emission from polarized ®Co. 
A backward-forward asymmetry was observed, in agreement with 
the results of Wu et al. (1957), thus confirming that parity is not conserved 
in B-decay. Although this detection method does not enable energy 
discrimination to be made, it has the merit of simplicity and of allowing 
the angular distribution to be measured simultaneously at different 
angles. In the present experiment, observations were confined to the 
0° and 180° directions. 

The photographic film has to be used within the low temperature 
cryostat at 1°K and therefore the main requirement of the film is that it 
should be sensitive to B-rays at this temperature. Several films were 
tested and most of these showed a large reduction in sensitivity as the 
temperature was lowered. However, it was found that for [ford Indus- 
trial G x-ray film, the blackening sensitivity to B-radiation from Co 


oe 
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is only reduced by a factor of two when the film is cooled from room 
temperature to 1°x. Most of the change in sensitivity occurs when 
the film is cooled from room temperature to 90°K. The y-rays have 
relatively little effect on the film owing to their much smaller absorption. 
Visible light was prevented from affecting the film by protecting it with 
aluminium foil 0-0002in. thick. Such a thickness reduces the Co 
B-ray intensity by about 10%. It was found that a two minute exposure 
at 1°K to a 10.0 source of Co at a distance of 5 mm produced sufficient 
blackening to detect a change greater than 4°/ in the intensity of the 
B-emission. (Scurlock and Sowter, to be published.) 


Section through demagnetization cryostat. A, radioactive crystal ; B, cylindri- 
cal photographic film ; C, apertures defining the B-ray beam ; D, alumi- 
nium window ; E, thread for raising the film carrier F ; G, spring guides 
for locating film carrier ; H, helium bath. 


The method of polarizing ®°Co was that described by Ambler ef al. 
(1953). A pair of single crystals of cerium magnesium nitrate was 
grown with a thin surface layer in which some of the Mg** ions were 
replaced by ®Co ions. The crystals were mounted in the demagnetization 
eryostat with the crystal axis horizontal. A sectional view of the appara- 
tus is given in the figure, and shows the cylindrical photographic film 
which could be raised past the crystals by an externally operated thread. 
The B-rays were screened from the film by a brass tube containing two 
apertures (4mm diameter) to define a fixed p-ray beam. 
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In the experiment, the crystals were cooled to 0-003°K by adiabatic 
demagnetization from 25 kilogauss applied perpendicularly to the crystal 
axis at 1°k, and the cobalt nuclei were polarized by applying a field 
of 600 gauss along the axis. The resulting temperature of the surface 
layer was of the order of 0-01°K and was estimated from the y-ray aniso- 
tropy measured with scintillation counters mounted outside the cryostat. 
The carrier was raised so as to expose the film for 2 minutes. Helium 
exchange gas was then admitted to warm the crystal to the temperature of 
the helium bath, and the carrier was raised again to take a normalizing 
exposure with the nuclei randomly oriented. Without exchange gas, 
the crystals would remain cold for nearly half an hour. The possibility 
of geometrical effects producing a 8-asymmetry was eliminated by repeat- 
ing the experiment with the polarizing field reversed. The apparatus 
was then warmed up to room temperature and the film removed and 
developed. The intensities of the exposures were measured with a 
microphotometer. 

Three sets of low temperature photographic runs with different pairs 
of crystals were made. All of these showed backward—forward asym- 
metries in which the greater intensity was against the direction of the 
polarizing field. The most consistent set showed an 18° difference 
between the 180° and 0° directions, with a corresponding y-ray anisotropy 
of 7%. The sign of the difference shows that the f-emission is preferred 
against the direction of the nuclear spin. 

If parity is not conserved, the angular distribution of the $-particles 
in an allowed transition may be expected to have the form 1+« cos @. 
(@ is the angle of emission measured from the axis of nuclear polarization.) 
The observed value of the asymmetry parameter « is —0-09+ 0-03 and the 
y-ray anisotropy (0-07) corresponds to a nuclear polarization P, of 0-25. 
The total error has been assessed from control experiments and arises 
from two main sources (1) a reduction in the B-sensitivity of the film by 
y-ray fogging, and (2) the difference in the experimental reproducibility 
of identical exposures. Since no energy discrimination is available, 
it is not possible to correct the experimental value of « for the effects 
of the magnetic field on the electron paths and of the back-scattering 
in the crystal. Our uncorrected result, however, agrees closely with that 
of Wu e¢ al. 

It is interesting to observe that the surface remains at about 0-01°K, 
as shown by the y-ray anisotropy, for a considerable time after demag- 
netization, and therefore the mechanism of thermal conduction is much 
more favourable than would be expected from an extrapolation of measure- 
ments on lattice conductivities made at higher temperatures. 
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Density Change of a Crystal Containing Dislocations 


By W. M. Lomser 


Atomic Energy Research Establishment, Harwell, near Didcot, Berks 
[Received June 5, 1957] 


SEEGER (1956) has recently evaluated the decrease in density of a crystal 
consequent upon the introduction of a dislocation. His calculation is 
based on the deviation of real crystals from ordinary linear elasticity. 
He estimates the magnitude of the effect by assuming that on shearing a 
crystal the planes of atoms ride up on the neighbouring planes to an 
extent proportional to the square of the shear strain, and determined 
the constant of proportionality analytically, using explicit interatomic 
forces. His result is that the expansion is about #2 per atomic plane of 
dislocation ({2=atomic volume), and that it is distributed in space just 
like the elastic energy (since both depend on the square of the elastic 
strains). 

It is interesting to note that the order of magnitude of the effect is 
confirmed by a very simple argument. The thermal expansion of most 
solids can also be related to the deviations from linear elasticity, and 
Gruneisen’s relation can be exhibited as a relation between thermal 
energy content and volume expansion. The Gruneisen constant y 1s 
defined as 

3aV 
gs OP 
where K is the. compressibility, « is the coefficient of linear expansion, 
and C,, is the specific heat. Rearranging and replacing « and U, by their 
usual definitions we get 


where H=thermal energy content. 
In a solid, half the thermal energy is kinetic energy, and half elastic 
potential energy. If the volume expansion is proportional to the elastic 


energy W, we finally get 
dV =2yKSW. 
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To substitute some numerical values, consider copper, with y=1-9, 
K=7 X10,” emt/ergs Vier 
bV 
— =2:6 10-¥ cm/er 
Sie oie [erg 
=4-25 x 10-*4 em?/ev. 


The energy of a dislocation in copper is typically about 6 ev/atom plane, 
and the atomic volume is 11-6 x 10-*4 cm®, so that we in fact predict a 
greater volume expansion than Seeger, nearly 22 per atomic length of 
dislocation line in copper. The discrepancy may be due to Seeger’s 
neglect of all but the shear strains parallel to the slip plane. 

It should be commented that the elastic constants effectively cancel in 
this calculation, since the energy of the dislocation is proportional to the 
elastic constants, and K inversely proportional. For an isotropic medium 
the volume expansion per unit length of edge dislocation line is 

(1—2¢) b? RB 


Yo 4r ae 
where o is Poisson’s Ratio, b the Burgers vector, R and ry the outer and 
inner radii of the effectively strained region. It would be interesting to 
compare experimentally the density recovery and energy release in say 
Cr (y=0-9) and Au and Pt (y=3-0). 
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REVIEWS OF BOOKS 


Atomic Physics. By Max Borw (6th edition revised in collaboration with 
R. J. Blin-Stoyle). Blackie & Son Ltd. [Pp. xiv+445.] Price 40s. 


Iv is twenty-two years since the first edition of this book appeared, and now 
the sixth comes to us looking as young and alluring as ever. The characteristic 
array of appendices remains, to satisfy appetites whetted by the equally 
characteristic and lucid text. There is a little more material than in previous 
editions, just enough to make us marvel at how little is needed to supplement 
what must have been, in the beginning, a most exciting and magisterial account 
of the new ideas in Physies. With familiarity some of the excitement has 
perhaps waned, but Born’s Atomic Physics still remains, and will long remain, 
we hope, the best introduction for the novice and the best corrective for the 
over-sophisticated initiate. JAN Se a 


Electricity and Magnetism. By B. I. BLnanny and B. Breaneny. (Clarendon 
Press.) [Pp. xiv+676.] Price 63s. 


Tuts excellent book by Prof. and Mrs. Bleaney supplies a real need in the 
advanced teaching of electricity and magnetism, for a book which shall be at 
once accurate, up-to-date and interesting (especially interesting). Some 450 
pages are devoted to those topics which may be expected to appear in such a 
book, and everywhere the account seems admirably intelligible, without undue 
emphasis on finer points of difficulty. There is even a photograph of a modern 
galvanometer, marking a notable break from traditional presentations. The 
last 200 pages treat phenomena in atomic and solid-state physics which are 
electric and magnetic in nature, such as fluctuations, dielectrics, metallic 
conduction, magnetism of solids, and magnetic resonance, and here the authors 
write with the authority of research workers in the field, while still bearing in 
mind the needs of the reader whom they are addressing. It is a pity that radio 
propagation and ionospheric physics are entirely neglected, but clearly this is 
not due to inadvertence. 

The book is full of examples and problems, many of which carry the argu- 
ments one stage further. The choice of examples is catholic and imaginative. 
One may confidently predict for it a long and useful life. ARB: 


Optics. By Bruno Rossi. (Addison-Wesley Publishing Co., Inc., Reading, 

Mass., U.S.A.) [Pp. 510.] Price $8.50. 

Tr is clear that, in this book, the author has given great thought to the best 
method of presentation of each topic. The style is clear and concise. This, 
together with the excellent diagrams and typography, makes the book a 
pleasure to read, and very easy to follow. 

The standard of the book is roughly that required for a normal degree course. 
It includes geometrical optics and physical optics. The chapters on the 
electromagnetic theory of light (which presuppose a knowledge of Maxwell's 
equations) would only be required by a student for final honours in Physics. 
For the rest of the book, the only mathematical knowledge which is assumed 
is a knowledge of calculus and simple differential equations. The student 
for final honours will find that the book does not cover a few special topics 
which he may require. These include the Fourier approach to diffraction, 
the theory of the resolving power of a microscope, and the theory of white 


light. 
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A large number of problems are given at the end of each chapter. These 
are more interesting than usual, and often extend the treatment given in the 
text. . 

For a first edition, the book is remarkably free from errors. One mistake, 
however, does occur on page 328, where it is stated that the transverse character 
of light waves cannot be explained on an elastic wave theory. 

This is a very valuable book, and will undoubtedly become a ‘standard 
textbook on the subject. Bate Ba 


The Journal of Photographic Science, March/Apru, 1957, Vol. 5, No. 2. (The 
Royal Photographic Society of Great Britain.) [Pp. xii+31.] Price: 
Annual subscription (six bi-monthly parts) £1 5s. Od. 


A sprcraL number of the Journal of Photographic Science has appeared dealing 
with high speed photography. The contents are as follows : 


“High frequency cinematography in the shock tube’ by H. Schardin ; 

“A 16mm. missile camera” by R. B. Herden ; 

‘* Photographic instrumentation at ordnance proving grounds” by G. 
Doughty ; 

‘ Detonation studies of pentolite sticks ’ by M. Sultanoff ; 

“High speed photography” by J. S. Courtney-Pratt and D. P. C. 
Thackeray ; 

‘Report of meeting : High-speed photography ” by P. Naslin. nore 

N. F. M. 


Les Dislocations. By J. Frrepei. (Gauthier-Villars, Paris.) [Pp. 314.] 

Price £3 14s. Od. 

Durine the last few years dislocation theory has undergone considerable 
development. The extent of this development can be judged from the scope 
of this book, the first part of which is concerned entirely with a purely formal 
treatment of dislocation theory. Chapters are included on elastic theory, 
glide and climb, on imperfect dislocations, and on crystal growth. The second 
part of the book is devoted to assemblies of dislocations and to applications 
of the theory to ‘ elastic’ moduli, work-hardening, recovery, recrystallization, 
creep and fracture. The last part of the book deals with the interaction of 
dislocations with impurities. 

The subject has been presented in a very clear manner ; often simple argu- 
ments are used to derive energies associated with dislocations and other defects. 
The order in which the theory is presented and the clear style make the book a 
valuable basis for a lecture course. On the other hand, numerous sections 
obviously represent new and original contributions to the subject. The 
section on climb is an example and should be read by anyone engaged in studies 
of high temperature creep and related phenomena. 

After reading this book one feels that whereas the basic properties of dis- 
locations are understood, much remains to be done in the field of application 
of the theory to the many-bodied problems of creep, work-hardening, re- 
crystallization etc. There is also a lack of calculations on the dislocation core, 
and this is emphasized by the fact that the Peierls-Nabarro treatment receives 
only cursory attention. 

This book is to be recommended to all those engaged in the study of dis- 
locations. It is by far the most extensive treatment of the subject yet produced. 
Unfortunately, circulation will however be restricted owing to the very high 
price. A cheaper translation into English would seem well worth while. 
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Physical Properties of Crystals. By J. F. Nye. (Clarendon Press, Oxford.) 

[Pp. 322+xv.] Price 50s. 

Tuts book deals with the systematic formulation of the macroscopic tensor 
properties of crystals. There are four parts: the first deals with the mathe- 
matics and properties of tensors. The second part covers the equilibrium 
properties of crystals, including magnetic susceptibility, polarization, piezo- 
electricity and elasticity, and ends with a chapter on the thermodynamic 
relations between these properties. Part III deals with transport properties, 
i.e. conductivity, and the last part with crystal optics. A number of useful 
appendices are added at the end. 

It has been a real pleasure to review this book, which has been written in a 
style of exemplary clarity. The arrangement of the subject matter is excellent 
and this makes the book, or rather parts of it, particularly useful as a basis 
for a lecture course. At the same time, most of the important tensor properties 
are treated in detail, and the book will therefore also serve as a useful reference 
book to solid state research workers. I believe that this book will remain 
the standard work on this subject for many years. Praise should also go to the 
Oxford University Press for the high quality of the production of this work. 
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Fig. 1 


A silver chloride crystal with surface parallel to the (001) planes after etching 
and exposure. Any dislocations which are present in these crystals 
cannot be made visible either by etching or by the internal separation of 
silver. The dark specks represent particles of silver on the surface 
of the crystal. (x 2500) 

Big. 2 


Dislocations forming elements of a coarse three dimensional network (see 
is) 

Hedges and Mitchell 1953, fig. 3) made visible in a crystal of silver 

chloride both by etching and by the internal separation of silver resulting 


This establishes the correspond- 


from the subsequent exposure to light. 
< 2500) 


ence between etch pits and dislocation lines. ( 


D. A. JONES and J. W. MITCHELL Phil. Mag. Ser. 8, Vol. 2, Pl. 40. 
Fig. 3 


A wall of dislocations in an annealed crystal of silver chloride. 
( x 2500) 


Dislocations introduced by the deformation of an element of the substructure of 
a crystal of silver chloride. (x 2500) 


Norts.—All the photomicrographs have been taken with phase contrast 


